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Introductory Remarks. 


The theoretical development* of the laws of motion of bodies is a problem of such interest 
and importance, that it has engaged the attention of all the most eminent mathematicians, 
since the invention of dynamics as a mathematical science by Galileo, and especially since the 
wonderful extension which was given to that science by Newton. Among the successors of 
those illustrious men, Lagrange has perhaps done more than any other analyst to give extent 
and harmony to such deductive researches, by showing that the most varied consequences 
respecting the motions of systems of bodies may be derived from one radical formula; the beauty 
ofthe method so suiting the dignity of the results, as to make of his great work a kind of scientific 
poem. But the science of force, or of power acting by law in space and time, has undergone 
already another revolution, ] and has become already more dynamic, by having almost dis- 
missed the conceptions of solidity and cohesion, and those other material ties, or geometrically 
imaginable conditions, which Lagrange so happily reasoned on, and by tending more and more 
to resolve all connexions and actions of bodies into attractions and repulsions of points: and 
while the science is advancing thus in one direction by the improvement of physical views, it 
may advance in another direction also by the invention of mathematical methods. And the 
method proposed in the present essay, for the deductive study of the motions of attracting or 
repelling systems, will perhaps be received with indulgence, as an attempt to assist in carrying 
forward so high an inquiry. 

In the methods commonly employed, the determination of the motion of a free point in 
space, under the influence of accelerating forces, depends on the integration of three equations 
in ordinary differentials of the second order; and the determination of the motions of a system 
of free points, attracting or repelling one another, depends on the integration of a system of 
such equations, in number threefold the number of the attracting or repelling points, unless we 
previously diminish by unity this latter number, by considering only relative motions. ‘Thus, in 
the solar system, when we consider only the mutual attractions of the sun and of the ten known 
planets,t the determination of the motions of the latter about the former is reduced, by the 
usual methods, to the integration of a system of thirty ordinary differential equations of the 
second order, between the coordinates and the time; or, by a transformation of Lagrange, $ to 
the integration of a system of sixty ordinary differential equations of the first order, between 
the time and the elliptic elements: by which integrations, the thirty varying coordinates, or the 
sixty varying elements, are to be found as functions of the time. In the method of the present 
essay, this problem is reduced to the search and differentiation of a single function, which 
satisfies two partial differential equations of the first order and of the second degree: and every 
other dynamical problem, respecting the motions of any system, however numerous, of attract- 

* [On the development of Theoretical Dynamics from Lagrange to Hamilton, see Cayley’s Report to the British 
Association, Dublin (1858), p. 1.] 

T [Hamilton is here referring to Boscovich's hypothesis on the constitution of matter (1743). Cf. Graves, Life 
of Hamilton, Vol. 1, p. 593, Vol. 1t, pp. 85-86, 88.] 

i [The planets known at the date of this paper were Mercury, Venus, Earth, Mars, Jupiter, Saturn, Uranus and 
four asteroids Ceres, Pallas, Vesta and Juno, the last being discovered in 1807. Neptune was not discovered until 


1846.] 
§ [Méc. Anal. 3rd ed. Tome rr, pp. 112-143.] 
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ing or repelling points, (even if we suppose those points restricted by any conditions of connexion 
consistent with the law of living force,) is reduced, in like manner, to the study of one central 
function, of which the form marks out and characterizes the properties of the moving system, 
and is to be determined by a pair cf partial differential equations of the first order, combined 
with some simple considerations. The difficulty is therefore at least transferred from the 
integration of many equations of one class to the integration of two of another: and even if it 
should be thought that no practical facility is gained, yet an intellectual pleasure may result 
from the reduction of the most complex and, probably, of all researches respecting the forces 
and motions of body, to the study of one characteristic function,* the unfolding of one central 
relation. 


The present essay does not pretend to treat fully of this extensive subject,—a task which 
may require the labours of many years and many minds; but only to suggest the thought and 
propose the path to others. Although, therefore, the method may be used in the most varied 
dynamical researches, it is at present only applied to the orbits and perturbations of a system 
with any laws of attraction or repulsion, and with one predominant mass or centre of pre- 
dominant energy; and only so far, even in this one research, as appears sufficient to make the 
principle itself understood. It may be mentioned here, that this dynamical principle is only 
another form of that idea which has already been applied to optics in the Theory of systems of 
rays, and that an intention of applying it to the motions of systems of bodies was announced} 
at the publication of that theory. And besides the idea itself, the manner of calculation also, 
which has been thus exemplified in the sciences of optics and dynamics, seems not confined to 
those two sciences, but capable of other applications; and the peculiar combination which it 
involves, of the principles of variations with those of partial differentials, for the determination 
and use of an important class of integrals, may constitute, when it shall be matured by the 
future labours of mathematicians, a separate branch of analysis.1 


Warrum R. HAMILTON. 
Observatory, Dublin, : 
March 1834. 


Integration of the Equations of Motion of a System, Characteristic Function of such 
Motion, and Law of varying Action. 

1. The known differential equations of motion of a system of free points, repelling or 
attracting one another according to any functions of their distances, and not disturbed by any 
foreign force, may be comprised in the following formula: 

X.m (x"ðx + y"dy 4-202) 2 9U. (1.) 
In this formula the sign of summation X extends to all the points of the system; m is, for any 


* Lagrange and, after him, Laplace and others, have employed a single function to express the different forces 
of a system, and so to form in an elegant manner the differential equations of its motion. By this conception, great 
simplicity has been given to the statement of the problem of dynamics; but the solution of that problem, or the 
expression of the motions themselves, and of their integrals, depends on a very different and hitherto unimagined 
function, as it is the purpose of this essay to show. 

T Transactions of the Royal Irish Academy, Vol. xv, p. 80. [Mathematical Papers, Vol. 1, p. 9.] A notice of 
this dynamical principle was also lately given in an article *On a general Method of expressing the Paths of 
Light and of the Planets," published in the Dublin University Review for October 1833. [Mathematical Papers, Vol. 
I, p. 311.] i 

į [Hamilton afterwards developed this suggestion in his Calculus of Principal Relations. See pp. 297—410.] 
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Secondly, the group, 
Bp» nito g ji Ms TR 
Sa, m05; $a, — Mzz; jas T, nin; 
òV EROR f òy y 
8b, = mbi; Na —n3b;; 3 — My by; (D.) 
ida AM ge gg iege ôy ;. 
8c, RM M €1; ca Moa Co; Scn — Manns 
and finally, the equation, 
oV 
sa (E.) 


So that if this function V were known, it would only remain to eliminate H between the 3» + 1 
equations (C.) and (E.), in order to obtain all the 3» intermediate integrals, or between (D.) 
and (E.) to obtain all the 3n final integrals of the differential equations of motion; that is, ulti- 
mately, to obtain the 3n sought relations between the 3» varying coordinates and the time, 
involving also the masses and the 6n initial data above mentioned; the discovery of which 
relations would be (as we have said) the general solution of the general problem of dynamics. 
We have, therefore, at least reduced that general problem to the search and differentiation of a 
single function V, which we shall call on this account the CHARACTERISTIC FUNCTION of motion 
of a system; and the equation (A.), expressing the fundamental law of its variation, we shall 
call the equation of the characteristic function, or the LAW OF VARYING ACTION. 


3. Toshow more clearly that the action or accumulated living force of a system, or in other 
words, the integral of the product of the living force by the element of the time, may be regarded 
as a function of the 6n + 1 quantities already mentioned, namely, of the initial and final co- 
ordinates, and of the quantity H, we may observe, that whatever depends on the manner and 
time of motion of the system may be considered as such a function; because the initial form of 
the law of living force, when combined with the 3» known or unknown relations between the 
time, the initial data, and the varying coordinates, will always furnish 3» + 1 relations, known 
or unknown, to connect the time and the initial components of velocities with the initial and 
final coordinates, and with H. Yet from not having formed the conception of the action as a 
function of this kind, the consequences that have been here deduced from the formula (A.) for 
the variation of that definite integral appear to have escaped the notice of Lagrange, and of the 
other illustrious analysts who have written on theoretical mechanics; although they were in 
possession of a formula for the variation of this integral not greatly differing from ours. For 
although Lagrange and others, in treating of the motion of a system, have shown that the 
variation of this definite integral vanishes when the extreme coordinates and the constant H 
are given, they appear to have deduced from this result only the well-known law of least action; 
namely, that if the points or bodies of a system be imagined to move from a given set of initial 
to a given set of final positions, not as they do nor even as they could move consistently with the 
general dynamical laws or differential equations of motion, but so as not to violate any supposed 
geometrical connexions, nor that one dynamical relation between velocities and configurations 
which constitutes the law of living force; and if, besides, this geometrically imaginable, but 
dynamically impossible motion, be made to differ infinitely little from the actual manner of 
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(z', y', 2’, being here, according to the analogy of our foregoing notation, the rectangular com- 
ponents of velocity of the point m, or the first differential coefficients of its coordinates taken 
with respect to the time;) an easy and well-known combination of the differential equations of 
motion, obtained by changing in the formula (1.) the variations to the differentials of the 
coordinates, may be expressed in the following manner, 


dT =dU, (5.) 
and gives, by integration, the celebrated law of living force, under the form 
T-U 4H. (6.) 


In this expression, which is one of the seven known integrals already mentioned, the quan- 
tity H is independent of the time, and does not alter in the passage of the points of the system 
from one set of positions to another. We have, for example, an initial equation of the same 
form, corresponding to the origin of time, which may be written thus, 

Ty —- U,4- H. (7.) 

The quantity H may, however, receive any arbitrary increment whatever, when we pass in 
thought from a system moving in one way, to the same system moving in another, with the 
same dynamical relations between the accelerations and positions ofits points, but with different 
initial data; but the increment of H, thus obtained, is evidently connected with the analogous 
increments of the functions 7' and U, by the relation 


AT — AU - AH, (8.) 
which, for the case of infinitesimal variations, may conveniently be written thus, 
oT —8U --8H; (9.) 


and this last relation, when multiplied by dt, and integrated, conducts to an important result. 
For it thus becomes, by (4.) and (1.), 


fz m (de Ba! dy By! de. Se) - [X .m (da: .d8a+dy’ . 99 -- dz 8z) + [SH dt, (10.) 


that is, by the principles of the calculus of variations, 


8V =X .m (x'8z + y'8y -- 2/02) — 21. m (a/8a + b'9b + c’5c) + 167, (A.) 
if we denote by V the integral 
t 
y- [x .m (x'da + y'dy + z'dz) =f 2T dt, (B.) 
0 


namely, the accumulated living force, called often the action of the system, frora its initial to 
its final position. 1 
If, then, we consider (as it is easy to see that we may) the action V as a function of the initial 
and final coordinates, and of the quantity H, we shall have, by (A.), the following groups of 
equations; first, the group, 


Ar SE pM sg BE 1 
x, 1*1» x, 2*"2? Òt, n"^n? 
NOn ins aibi io cleo idus oot (Q.) 
994 1413 8y, 23/2; ses y, nn: ^ 
èy j.p Ji èy 
8z, = 0421; 82, = iia 25; &, TIU, E, 
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Secondly, the group, 


Bp. sm a’: BPE '. Lini pee a’: 
Sa, Mı 1: $a, Moala; eee San T n "n? 
èy PNE T i oy I 
5b, — mbi; 5b, — Maba; + 8b. —m4b,; (D.) 
"s ci: iu PE ci: aga c!: 
óc Thee Se Meee iud BUH Ss NUI 
and finally, the equation, 
èy 
sgt (E.) 


So that if this function V were known, it would only remain to eliminate H between the 3n + 1 
equations (C.) and (E.), in order to obtain all the 3n intermediate integrals, or between (D.) 
and (E.) to obtain all the 3n final integrals of the differential equations of motion; that is, ulti- 
mately, to obtain the 3» sought relations between the 3» varying coordinates and the time, 
involving also the masses and the 6n initial data above mentioned; the discovery of which 
relations would be (as we have said) the general solution of the general problem of dynamics. 
We have, therefore, at least reduced that general problem to the search and differentiation of a 
single function V, which we shall call on this account the CHARACTERISTIC FUNCTION of motion 
of a system; and the equation (A.), expressing the fundamental law of its variation, we shall 
call the equation of the characteristic function, or the LAW OF VARYING ACTION. 


3. Toshow more clearly that the action or accumulated living force of a system, or in other 
words, the integral of the product of the living force by the element of the time, may be regarded 
as a function of the 6n+1 quantities already mentioned, namely, of the initial and final co- 
ordinates, and of the quantity H, we may observe, that whatever depends on the manner and 
time of motion of the system may be considered as such a function; because the initial form of 
the law of living force, when combined with the 3n known or unknown relations between the 
time, the initial data, and the varying coordinates, will always furnish 3» 4- 1 relations, known 
or unknown, to connect the time and the initial components of velocities with the initial and 
final coordinates, and with H. Yet from not having formed the conception of the action as a 
function of this kind, the consequences that have been here deduced from the formula (A.) for 
the variation of that definite integral appear to have escaped the notice of Lagrange, and of the 
other illustrious analysts who have written on theoretical mechanics; although they were in 
possession of a formula for the variation of this integral not greatly differing from ours. For 
although Lagrange and others, in treating of the motion of a system, have shown that the 
variation of this definite integral vanishes when the extreme coordinates and the constant H 
are given, they appear to have deduced from this result only the well-known law of least action; 
namely, that if the points or bodies of a system be imagined to move from a given set of initial 
to a given set of final positions, not as they do nor even as they could move consistently with the 
general dynamical laws or differential equations of motion, but so as not to violate any supposed 
geometrical connexions, nor that one dynamical relation between velocities and configurations 
which constitutes the law of living force; and if, besides, this geometrically imaginable, but 
dynamically impossible motion, be made to differ infinitely little from the actual manner of 
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motion of the system, between the given extreme positions; then the varied value of the 
definite integral called action, or the accumulated living force of the system in the motion thus 
imagined, will differ infinitely less from the actual value of that integral. But when this well- 
known law of least, or as it might be better called, of stationary action, is applied to the deter- 
mination of the actual motion of a system, it serves only to form, by the rules of the calculus of 
variations, the differential equations of motion of the second order, which can always be 
otherwise found. It seems, therefore, to be with reason that Lagrange*, Laplace, and Poisson 
have spoken lightly of the utility of this principle in the present state of dynamics. A different 
estimate, perhaps, will be formed of that other principle which has been introduced in the present 
paper, under the name of the law of varying action, in which we pass from an actual motion to 
another motion dynamically possible, by varying the extreme positions of the system, and (in 
general) the quantity H, and which serves to express, by means of a single function, not the 
mere differential equations of motion, but their intermediate and their final integrals. 


Verifications of the foregoing Integrals. 


4. A verification, which ought not to be neglected, and at the same time an illustration of 
this new principle, may be obtained by deducing the known differential equations of motion 
from our system of intermediate integrals, and by showing the consistence of these again with 
our final integral system. As preliminary to such verification, it is useful to observe that the 
final equation (6.) of living force, when combined with the system (C.), takes this new form, f 


ON " 


and that the initial equation (7.) of living force becomes by (D.) 


STRONE e 


These two partial differential equations, initial and final, of the first order and the second 
degree, must both be identically satisfied by the characteristic function V: they furnish (as we 
shall find) the principal means of discovering the form of that function, and are of essential 
importance in its theory.t If the form of this function were known, we might eliminate 3n — 1 
of the 3n initial coordinates between the 3n equations (C.); and although we cannot yet perform 
`- the actual process of this elimination, we are entitled to assert that it would remove along with 
the others the remaining initial coordinate, and would conduct to the equation (6.) of final 
living force, which might then be transformed into the equation (F.). In like manner we may 
conclude that all the 3n final coordinates could be eliminated together from the 3» equations 
(D.), and that the result would be the initial equation (7.) of living force, or the transformed 
equation (G.). We may therefore consider the law of living force, which assisted us in discovering 
the properties of our characteristic function V, as included reciprocally in those properties, and 
as resulting by elimination, in every particular case, from the systems (C.) and (D.); and in 
treating of either of these systems, or in conducting any other dynamical investigation by the 


* [For Lagrange's remarks on the principle of least action, see Méc. Anal. 3rd ed. Tome 1, pp. 229, 230.] 

T [These equations are exactly analogous to the equations Q=0, Q’=0 of the Third Supplement. Mathematical 
Papers, Vol. 1, pp. 170, 485.] 

i [For Jacobi's criticism, see Appendix, p. 613.] 
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method of this characteristic function, we are at liberty to employ the partial differential 
equations (F.) and (G.), which that function must necessarily satisfy. 


It will now be easy to deduce, as we proposed, the known equations of motion (3.) of the 
second order, by differentiation and elimination of constants, from our intermediate integral 
system (C.), (E.), or even from a part of that system, namely, from the group (C.), when com- 
bined with the equation (F.). For we thus obtain 
doy. ov... oF. NW. 
die 3 Sag 75546 7 Ee. 

ôy Lá ; 21V 
HH Sman Y may +4" Sandy, 
+2, èy AT. ôy j èy 
82,52, |? da "wn ^ $2, 02, 

ló5V eV 15V Py l3, WF 
Om, Bx, Bi tm, Se, Serd, 7 m, Bx, Bx, Sz, io 
LUN ME ASEE Aa 

m,9y,0z,8y, m48y,0z,8y, ^" m, SY, 52, dy, 

15V &V 18V Py 1 8V OM 


r ipw RM er T Sz, 52,52, 


$^ Zw ic) +h) *()]- = 55, (UH) 


m,2z;- = ^ (12.) 
And in like manner we might deduce, by differentiation, from the integrals (C.) and from (F.) 
all the other known differential equations of motion, of the second order, contained in the set 
marked (3.); or, more concisely, we may deduce at once the formula (1.), which contains all 
those known equations, by observing that the intermediate integrals (C.), when combined with 
the relation (F.), give* 


m,xi= 


that is, we obtain 


TAUN AE: ela à id z 
Vè 8V8 VÒ oV az 
"Pu m Su Sat By 3g Ba Sa) A a =”) 


=E (te 55+ 805, +852) -am Css) + (59) * (ae) | ue 


-X(wie5miemr 5) (0 2 
-8U. 


5. Again, we were to show that our intermediate integral system, composed of the equations | 
(C.) and (E.), with the 3» arbitrary constants a, , b1, c,, ... Ans Ons Cn, (and involving also the 


* [rette i Sm E ioi APES i ae pe, Ee ms (iz) | Ome P ; b= 80, | 
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auxiliary constant H,) is consistent with our final integral system of equations (D.) and (E.), 
which contain 3n other arbitrary constants, namely, aj, bj, cj, ... a), bi, cp. The immediate 
differentials of the equations (C.), (D.), (E.), taken with respect to the time, are, for the first 


group, 
eee ea d àV 


dism, VE dt8n, ee di&z, nn; 
d oV NW i14 ii doV | ý 
di dy, ai diBy, aye; T diBy, nn: (H.) 
LE WWW opp A uL TT 
di, v gig, Cv 7 di, e 
for the second group, 
EET UD TIT et Lb: as 
QM o NS 160r gv 
doV 1-4 P 4 GOV). -a 
dio, e. dt db, 0; eee a, °° (I.) 
ALLEE NR 1.1 AN T 
EN ENS s oM. 
and finally, for the last equation, 
dV 
dsH b (K.) 


By combining the equations (C.) with their differentials (H.), and with the relation (F.), we 
deduced, in the foregoing number, the known equations of motion (3.); and we are now to show 
the consistence of the same intermediate integrals (C.) with the group of differentials (I.), 
which have been deduced from the final integrals. 


The first equation of the group (I.) may be developed thus: 


oaa Ei a Eis ua n 
1$9a,02,..*0a,0m, °°." Sa, 8a, 
Pe) xo dii , &y 
tY Sa Sy, ay, "^57" 3a, dy, (14.) 
PAT op CABE A iao fis) oM Eo 
[Mae TOES ui 000," 


and the others may be similarly developed. In order, therefore, to show that they are satisfied 
by the group (C.), it is sufficient to prove that the following equations are true, 


RR 

o err Cry E : 
2 2 2 

t enr or CE) | 


the integer 7 receiving any value from 1 to n inclusive; which may be shown at once, and the 
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required verification thereby be obtained, if we merely take the variation of the relation (F.) 
with respect to the initial coordinates, as in the former verification we took its variation with 
respect to the final coordinates, and so obtained results which agreed with the known equations 
of motion, and which may be thus collected, 


Pian qu MEAE 
ôx, 2m |\ 8a y oz) ) Sx, 
8 l((8VY (dV\?_ [80V 8U 
s 2m (ae) (59) e) Jo 
2g Lf e 
àz; sats ap Ate) joa 

The same relation (F.), by being varied with respect to the quantity H, conducts to the 


expression 
ò REE E A E E A E i A 
se? an ss *(s;) (5) |=1; TA 


and this, when developed, agrees with the equation (K.), which is a new verification of the 
consistence of our foregoing results. Nor would it have been much more difficult, by the help 
of the foregoing principles, to have integrated directly our integrals of the first order, and so to 
have deduced in a different way our final integral system. 


6. It may be considered as still another verification of our own general integral equations, 
to show that they include not only the known law of living force, or the integral expressing that 
law, but also the six other known integrals of the first order, which contain the law of motion 
of the centre of gravity, and the law of description of areas. For this purpose, it is only necessary 
to observe that it evidently follows from the conception of our characteristic function V, that 
this function depends on the initial and final positions of the attracting or repelling points of a 
system, not as referred to any foreign standard, but only as compared with one another; and 
therefore that this function will not vary, if without making any real change in either initial or 
final configuration, or in the relation of these to each other, we alter at once all the initial and 
all the final positions of the points of the system, by any common motion, whether of translation 
or of rotation.* Now by considering three coordinate translations, we obtain the three following 
partial differential equations of the first order, which the function V must satisfy, 


oV oV 

3173.35 
oV aV 
oV oV 

y ial cabe 


and by considering three coordinate rotations, we obtain these three other relations between 
the partial differential coefficients of the same order of the same characteristic function, 


* ['The function V is obviously independent of the choice of coordinate axes and the motions considered are 
equivalent to changes in the coordinate axes. Hence the value of V does not alter under such motions.] 
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eV OV SF. 38 


oV eV oV oV 
ST" car 8E a BO 
zr) 23s 5) —° 


and if we change the final [differential] coefficients of V to the final components of momentum, 
and the initial coefficients to the initial components taken negatively, according to the 
dynamical properties of this function expressed by the integrals (C.) and (D.), we shall change 
these partial differential equations (O.), (P.), to the following, 


àE.mz'-ià.ma'; ZX.my'2Z.mb'; X.mz'-2Z.m«c; (15.) 

and à .m (xy’ —yx') 2 XA.m (ab' —ba’); 
X.m(yz —2y'):X.m oc -an;| (16.) 

X .m (2x’ — x2’) 2 XA. m (ca/ — ac’). 


In this manner, therefore, we can deduce from the properties of our characteristic function 
the six other known integrals above mentioned, in addition to that seventh which contains the 
law of living force, and which assisted in the discovery of our method. 


Introduction of relative or polar Coordinates, or other marks of position of a System. 


7. The property of our characteristic function, by which it depends only on the internal or 
mutual relations between the positions initial and final of the points of an attracting or repelling 
system, suggests an advantage in employing internal or relative coordinates; and from the 
analogy of other applications of algebraical methods to researches of a geometrical kind, it may 
be expected that polar and other marks of position will also often be found useful. Supposing, 
therefore, that the 3n final coordinates z, , Y1, 2; , ... Ens Yn» Zn have been expressed as functions 
of 3n other variables, 5, , 72; ... Nan, and that the 3» initial coordinates have in like manner been 
expressed as functions of 3n similar quantities, which we shall call e; , e;, ... egn, we shall proceed 
to assign a general method for introducing these new marks of position into the expressions of 
our fundamental relations. 

For this purpose we have only to transform the law of varying action, or the fundamental 
formula (A.), by transforming the two sums, 

X.m(a'8rz--y/9y--202), and  X.m(a'8a 4 b'b - c'8c), 
which it involves, and which are respectively equivalent to the following more developed 
expressions, 
E.m (82: 4- y'8y - 202) 2 m, (x, 8a, +y; Sy, +2182) 
+ Mz (22 3%_ t y39ya -- 25925) (17.) 
 &c. + m, (Ep 9v, t Yn 9y,, + 2,92,); 
X. m (a'a + b'8b + c'8c) = m, (a, 8a, 4- 559b, 4-190) 
+ m (a5 da, +63 8b, + c58c,) (18.) 
T &c. + m, (a, 9a, + b; db, + Cn 8c, ). 
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Now 2; being by supposition a function of the 3n new marks of position 7, , ... 137? He vacation 
82, , and its differential coefficient z;, may be thus expressed: 
ox ba 19. 
aid 7 hey na +++ AM vs 9035; ( 
, àv; , da; , A , (20.) 
X Sn, tae Net. Yi Nn» 


and similarly for y; and z;. If, then, we consider x; as a function, by (20.), of, .-- Nn» involving 
also in general 7, , ... nsn» and if we take its partial differential coefficients of the first order with 


respect to 71, ... Ngn, We find the relations, 
Sar, Ox, 8x; Om, Or; Ox, (21.) 
Sn ony Sna ông’ ie 895, E 


and therefore we obtain these new expressions for the variations 9z;, dy;, 52;, 


òx; òx; òx; 
_ by w i 
d ym x UL +e UT let. a nsn» (22.) 


D i Òz; 
da= 51 h tp at eee gym 


Substituting these expressions (22.) for the variations in the sum (17.), we easily transform 
it into the following, 


E.m (z'9z - y'8y - 2/82) 2 X. m(2’ Dail -+ y» iz) .O, 
ôni 1 èni 


àz' 8y 82’ 
+E.m(2 Ty! <5 zx ^ 
05 , 5 873 9n. 


òx’ y' 82’ 
+ &e.4E.m(2’ ; d &z)- 
"E Y y. Ton 


(23.) 


zx $i S IW deti. ET 903,; 


T being the same quantity as before, namely, the half of the final living force of the system, but 
being now considered as a function of 7}, ... Ngn, involving also the masses, and in general 
7» +++ Ngn» and obtained by substituting for the quantities a’, y’, z’ their values of the form (20.) 
in the equation of definition 
T =42.m(a'2+y'24+2). (4.) 
In like manner we find this transformation for the sum (18.), 


X. m (õa +b o'e) = e E + 52 Bea. (24.) 
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The law of varying action, or the formula (A.), becomes therefore, when expressed by the 


present more general coordinates or marks of position, * 
BV -X. A7 x. da Be 18H; (Q.) 


and instead of the groups (C.) and (D.), into which, along with the equation (E.), this law 
resolved itself before, it gives now these other groups, 


or. OT OP ST. I LA. (R.) 
905, «y Sy, Sna "^ Snan «y 
and 
M. HK Mss ELM, WS 
NS X Uu vod MS Uu. ; 


The quantities e, , €z, ... €gn, and €j, €z, ... €5,,, are now the initial data respecting the manner 
of motion of the system; and the 3n final integrals, connecting these 6n initial data, and the 
n masses, with the time ż, and with the 3n final or varying quantities 74, 72, ... Ngn, Which mark 
the varying positions of the » moving points of the system, are now to be obtained by eliminating 
the auxiliary constant H between the 3n + 1 equations (S.) and (E.); while the 3» intermediate 
integrals, or integrals of the first order, which connect the same varying marks of position and 
their first differential coefficients with the time, the masses, and the initial marks of position, 
are the result of elimination of the same auxiliary constant H between the equations (R.) and 
(E.). Our fundamental formula, and intermediate and final integrals, can therefore be very 
simply expressed with any new sets of coordinates; and the partial differential equations (F.), 
(G.), which our characteristic function V must satisfy, and which are, as we have said, essentia! 
in the theory of that function, can aiso easily be expressed with any such transformed co- 
ordinates, by merely combining the final and initiai expressions of the law of living force, 


T=U+H, (6.) 

T,=U,+H, (7.) 
with the new groups (R.) and (S.). For this purpose we must now consider the function U, of 
the masses and mutual distances of the several points of the system, as depending on the new 


* [Equation (Q.) can be obtained most simply as follows without introducing the cartesian coordinates. 
Let T be a function of the ys and 7s, being homogeneous of the second degree in the latter, and satisfying 
T — U +H, where U is a function of the n’s, and H a constant along each of the curves considered. Writing 


tar 
v=2 f’ riaj ba 
ay 


we have èy 2 (5) f EQ. ds 
ar ar (ar 
ae ste mB) anf uC) 
anf jj) =j, see) ™ 
But ' Ta =2T=T+U+H, 
ar 
and therefore zT JL ‘=X inte +8H, 
By) cpm 
im, dàT àT bU 
so that we find BV =E pon- E G be te [ET 5 5) onde 


This, by the principle of least action, leads to the equations (Z.) and equation (Q.) follows immediately.] 
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marks of position 7,, 72, ... Ngn; and the analogous function Ug, as depending similarly on the 
initial quantities e, , €g, ... €3n; we must also suppose that T'is expressed (as it may) as a function 


of its own coefficients i ; i we A. , which will always be, with respect to these, homogeneous 
i 2 3n 
of the second dimension, and may also involve explicitly the quantities 71, N2, ... Nan; and that 
To is expressed as a similar function of its coefficients a ais. HE oT». so that 
2 n 
1 2 7]3n (25.) 
T,-F (sz $e nne ) 
Met Bey” Beha? 


and that then these coefficients of T' and T', are changed to their values (R.) and (S.), so as to 
give, instead of (F.) and (G.), two other transformed equations, namely, 


eV èy an 
wT) wo bee waca =U+H, T. 
Bm s 
and, on account of the homogeneity and dimension of To, 
eV àV . 8V 
P (s fegh je) - Ue H. (U.) 


8. Nor is there any difficulty in deducing analogous transformations for the known differ- 
ential equations of motion of the second order, of any system of free points, by taking the 
variation of the new form (T.) of the law of living force, and by attending to the dynamical 
meanings of the coefficients of our characteristic function. For if we observe that the final 
living force 27', when considered as a function of 7, na, ... Ngn» and of 51, 3, ... Ngn» IS neces- 
sarily homogeneous of the second dimension with respect to the latter set of variables, and must 
therefore satisfy the condition 


eT . ,8T oT 
2T =ni s sedan ` 26. 
7A i N2 Sn Ur Stan (26.) 
we shall perceive that its total variation, 
oT oT 
òT — — 0n t ;—99, VE RE 
90, dn» 753;, (27.) 
Th A cl a + POR Sni | 
ur Moni; Sie eoe 
may be put under the form 
oT oT oT 
óT =718— +7999 <5 +e 9349 —;- 
7 wd Ur wa Nan + 
oT oT oT 
—<— 05, — —— 894 — Snan 
8n na 973, 
ST 8 GP) 
=> a 
185 5 71 
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and therefore that the total variation of the new partial differential equation (T.) may be thus 
written, 


(oe òT U 


in which, if we observe that 7’ = y and that the quantities of the form » are the only ones 
which vary with the time, we shall see that 


àV (d 8V d èV d èV 
Bait Ge -X[as. 1+ Se be) + + 35H 


because the identical equation 8dV —d98V gives, when developed, 
GAA apa: deor. dH 


.6H, (29.) 


ò 8e oH ` 


V V 
LET be) edis. SH. 


Decomposing, therefore, the expression (V.), for the variation of half the living force, into as 
many separate equations as it contains independent variations, we obtain, not only the equation 


-z(a T (30.) 


d 8V 
which had already presented itself, and the group 
d 8V d èV — d èV 
Bu gaT h T AW) 


which might have been at once obtained by Miu ee from the final integrals (S.), but also 
a group of 3» other equations of the form 


dóV 8T U 
diy 5 ~ òn T 
which give, by the intermediate integrals (R.), 
PIE 87 8U. 
diby m By NER 


that is, more fully, 
dar BT 3U 


dt Sna Pw, bey (2) 


...... 


These last transformations of the differential equations of motion of the second order, of an 
attracting or repelling system, coincide in all respects (a slight difference of notation excepted, ) 
with the elegant canonical forms in the Mécanique Analytique of Lagrange*; but it seemed 


* [Méc. Anal. 3rd ed. Tome 1, pp. 290-292.] 
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worth while to deduce them here anew, from the properties of our characteristic function. And 
if we were to suppose (as it has often been thought convenient and even necessary to do,) that 
the » points of a system are not entirely free, nor subject only to their own mutual attractions 
or repulsions, but connected by any geometrical conditions, and influenced by any foreign 
agencies, consistent with the law of conservation of living force; so that the number of indepen- 
dent marks of position should be now less numerous, and the force-function U less simple than 
before; it might still be proved, by a reasoning very similar to the foregoing, that on these sup- 
positions also (which, however, the dynamical spirit is tending more and more to exclude,*) 
the accumulated living force or action V of the system is a characteristic motion-function of the 
kind already explained; having the same law and formula of variation, which are susceptible of 
the same transformations; obliged to satisfy in the same way a final and an initial relation be- 
tween its partial differential coefficients of the first order; conducting, by the variation of one of 
these two relations, to the same canonical forms assigned by Lagrange for the differential 
equations of motion; and furnishing, on the same principles as before, their intermediate and 
their final integrals. To those imaginable cases, indeed, in which the law of living force no longer 
holds, our method also would not apply; but it appears to be the growing conviction of the 
persons who have meditated the most profoundly on the mathematical dynamics of the 
universe, that these are cases suggested by insufficient views of the mutual actions of body. 


9. It results from the foregoing remarks, that in order to apply our method of the cha- 
racteristic function to any problem of dynamics respecting any moving system, the known law 
of living force is to be combined with our law of varying action; and that the general expression 
of this latter law is to be obtained in the following manner. We are first to express the quantity T', 
namely, the half of the living force of the system, as a function (which will always be homo- 
geneous of the second dimension,) of the differential coefficients or rates of increase 7; , 75, &c., 
of any rectangular coordinates, or other marks of position of the system: we are next to take the 
variation of this homogeneous function with respect to those rates of increase, and to change 
the variations of those rates 97; , 955, &c., to the variations 87, , ôn, &c., of the marks of position 
themselves; and then to subtract the initial from the final value of the result, and to equate the 
remainder to èV —18H.T A slight consideration will show that this general rule or process for 
obtaining the variation of the characteristic function V is applicable even when the marks of 
position 71, 72, &c. are not all independent of each other; which will happen when they have 
been made, from any motive of convenience, more numerous than the rectangular coordinates 
of the several points of the system. For if we suppose that the 3» rectangular coordinates 
24, Y1 Zis + Uns Yn» Zn have been expressed by any transformation as functions of 3n + k other 
marks of position, 14, 72, «+» 93444; Which must therefore be connected by k equations of 
condition, 

O= 4 (01.912; +++ 192; 
0 — ds (m1 N2 +++ Nan+k)> (31.) 


02d. (m 75; --- Nsn+k)> 


* [See Reference to Boscovich, p. 104.] 
t Ee By -E e= 7 ui. | 
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giving k of the new marks of position as functions of the remaining 3n, 


N3n+1 = Ja (m T2577 Nan)» 
7)3n4-2 — be (71572, --- Nan)» 


(32.) 
7)3n--k = Pr (3 N2» ++ Nan)» 
the expression T =42.m(x'2+y'2 +z?) (4.) 


will become, by the introduction of these new variables, a homogeneous function of the second 
dimension of the 3n + k rates of increase 7} , 73, ... 95,4, involving also in general 7, , 7, 


*t*7]3n--k 
and having a variation which may be thus expressed: 
oT oT oT ; 
oT = (= ; oni + (s ; ông tec [or -) ÒNsn+k 
tf (5) an+ 5 (Z) Siti: + oT IL i (33. 
i ? òn die ón ^ O75, ik ie ipi 
- or in this other way, 
oT oT oT 
oT = M az; Oi + ôn By omnt ‘Ding. $ Onan 
E. ST 57" 68.) 


rt Se PETE. 09g t .. Eu ~ Stan» J 


on account of the relations (32.), Das tr when Re with respect to the time, 


à 
Nsn+1= nf 9h » i Ne n Tec "s oot 89g Pa 
M. ; 9s : M. 
Nen+2= "Ux 80, +73 89, Tec UMP ee 875; (35.) 
Op, , by [ON 


Nsn+k = ees Sn LT Fee FY Loo 
pen therefore, attending only to the variations of quantities of the form 7’, 


à 
9054,44 = st nit Bet at Dg 755 ; 
’ ò 8 
E Ha m+ e sg 975, ; (36.) 
kayode 5 x ^ | 
On, Li = s oni T 2 òns Ne pore 85; Ye ha 


Comparing the two ota (33.) and ies we find is (36.) the relations 


eT - (55 (i | oU, s | y, ( eT | ob, 
"d NEM T TC 7 [ow 
õn; 0m; ‘ ji 804 eig 8m 9735.Lk ôn 


AEO 
855 003 903441 93 0035.12 O1 2 SN3n+k 13 (37.) 


oT 2 ( òT E oT ) oT IE: 
aT Dor QUI IC ERN +l<7 zÉ et IRL Ero 
n5; (Se 993541 Ngn 9035.12 nsn men 90s, 


www.rcin.org.pl 


120 II. GENERAL METHOD IN DYNAMICS [9 
which give, by (32.), 


iom 57 Ona. EL 7 7- Šan = (=) 09, + =) na +.. oam Jon (38.) 


we may therefore put the expression (Q.) under the following more general form, 


pr ze) 8e +tôH, (A1) 


v -z.( 


the coefficients (= 3) being formed by treating all the 3n +k quantities 71, 75, ... nin. aS 


òn’ 
independent; which was the extension above announced, of the rule for forming the variatiòn 
of the characteristic function V. 

We cannot, however, immediately decompose this new expression (A!.) for 8V, as we did 
the expression (Q.), by treating all the variations 85, 9e as independent; but we may decompose 
it so, if we previously combine it with the final equations of condition (31.), and with the 
analogous initial equations of condition, namely, 


0=9, (êi êz, eee Can+k)> 


0= o, (e ;€95 eee bend) (39.) 


0= D; (e1, Cases Can+k)» 


which we may do by adding the variations of the connecting functions ¢,, ... $4, Py, ... Bp, 
multiplied respectively by factors to be determined, A,, ... Aj, A,, ... Ap. In this manner the 
law of varying action takes this new form, 


8y-E. (i) - z. (52) seam ex. ABj E. A80; (BL) 


and decomposes itself into 6n + 2k + 1 separate expressions, for the partial differential coeffi- 
cients of the first order of the characteristic function V, namely, into the following, 


oV -( 2) OP x a X. 96s 89, 
+A +...4A,——, 
Sn om 4 8 aht. a 1 òn 
SV _ (82), y ipa, Ms A 
57; * 13g, Pago Fog (C1) 
èy ( òT | 89, E 
=(<7—] +à +... +r; ; 
ÒNan+k ÒNsn+k 1 Onsn ik kon 
and 
oV — P oT, 806, 96, 89, ‘ 
8e (Set) ettet 
oV — 9T, 86, 90, 89, 
8e, (s) + Avge iwi (D1.) 
8V ( 8T, | 86, 8d, 
= -| ] +45 +. +A : 
Oan. Lk O65 nik y Seansk l s 6s, ik 
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besides the old equation (E.). The analogous introduction of multipliers in the canonical forms 
of Lagrange, for the differential equations of motion of the second order, by which a sum such as 


x Ag is added to in the second member of the formula (Y.), is also easily justified on the 


principles of the present essay. 


Separation of the relative motion of a system from the motion of its centre of gravity ; 
characteristic function for such relative motion, and law of its variation. 


10. As an example of the foregoing transformations, and at the same time as an important 
application, we shall now introduce relative coordinates, x, , y, , z, , referred to an internal origin 
£,» Y, 2,; that is, we shall put 


CiS ta tEn, YeqVYutyn, BAZ t2,» (40.) 
and in like manner 
b,=b,,+6 (41.) 


Q4 —0,,T4,, n? C=C, tC; 


together with the differentiated expressions 

XP—X44T, =Y tY, TJ -—*]4tt, (42.) 
and 

a,=a),ta,, b,=b),+b', C=C ter. (43.) 


"? n? 


Introducing the expressions (42.) for the rectangular components of velocity, we find that the 
value given by (4.) for the living force 27! decomposes itself into the three following parts, 


2T =X .m (x2 + y'2+2'2) Z Yi.m(z?--y,?-2;?) 
T2 (x, X.mz,-ty,2X.my,T2,2.mz,)-4 (2,?- y;? 4-2,?) Xm; (44.) 


if then we establish, as we may, the three equations of condition, 


x.mz,=0, L.my,=0, 2ZX.m2,-0, (45.) 
which give by (40.), 
à.mz x.my x.mz 
Tam Sg? es. om kt Ems. at 


so that xz, , y,, z, are now the coordinates of the point which is called the centre of gravity of 
the system, we may reduce the function 7' to the form 


TaT-T. (47.) 
in which 
T, 24X.m(v? ry? rz?) (48.) 
and 
T,=} (2; +Y, +2,) Um. (49.) 


By this known decomposition, the whole living force 2T of the system is resolved into the 
two parts 27', and 27, , of which the former, 27',, may be called the relative living force, being 
that which results solely from the relative velocities of the points of the system, in their motions 
about their common centre of gravity x, , y, , z,; while the latter part, 27', , results only from the 
absolute motion of that centre of gravity in space, and is the same as if all the masses of the 
system were united in that common centre. At the same time, the law of living force, 
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T — U -- H, (6.), resolves itself by the law of motion of the centre of gravity into the two 
following separate equations, | 


T,=U+H,, (50.) 
and 
T  =H,; (51.) 
H, and H, being two new constants independent of the time ż, and such that their sum 
H,+H,=H. (52.) 


And we may in like manner decompose the action, or accumulated living force V, which is 


t 
equal to the definite integral [ 2T dt, into the two following analogous parts, 
Jo 


V=V,+7,, (E) 
determined by the two equations, 
pi Í ' 2T. di, (FL) 
and ^ 
Vial | ‘oT dt. (G1) 
The last equation gives by (51.), t i 
V, 22H,t; (53.) 
a result which, by the law of motion of the centre of gravity, may be thus expressed, * 
V, =V (x, —a,)* + (y, —6,) + (2, —c,* . V 2H, Em: (Hi.) 
a, , b,, c, being the initial coordinates of the centre of gravity, so that 
aja Eee, baie, o EP (54.) 


And for the variation 5V of the whole function V, the rule of the last number gives 
ôV =E .m (x dx,—a’ da, +y’ dy, — b, db, - 2,82, —c,0c,) 
+(x’ a, — a; da, -- y, dy, — b; 8b, +2) 82, — c; dc,) Lm 
+t8H +A, X.msx,+A,X.msy, AX mz, 
+A,X.mda,+A,&.mdb,+A,2X.méc,; 


while the variation of the part V,, determined by the equation (H!.), is easily shown to be 
equivalent to the part 


òV, = (x; x, — a, da, +Y, dy, — b, 5b, -- 2, dz, — c, 9c, ) m + tH; (K1.) 
the variation of the other part V, may therefore be thus expressed, 
SV, 2 X.m (x du, —a;9a, 4- y, dy, — b; òb, +z, dz, — c, 9c,) 
+t8H, +A, X.msx,+A,X.mby, +AsX.msz, (L.) 
TAÀAQàE.mà$O08,-- A,€3.m0b,-- A,£2.m8e,: 
and it resolves itself into the following separate expressions, in which the part V, is considered 


(I*.) 


* [If v is the velocity of the centre of gravity, 
(x, —a,)* + (Yn — Oy)? + (Zn — Co)? v 2 2H, |m.] 
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as a function of the 6n + 1 quantities z,;, y,;, 2,;, 9,;, 6,;,¢,;, H, , of which, however, only 6n — 5 


are really independent: 


first group, 
T =M HAM; .. af =Mn Xin, xm 
za 02, 
P. = my yi £m; S maya Ama; 
8y4 12;1 s EL 8y,, n9 n 2 n? 
Pihaus : 2L, a. Zin T Ag; 
8g. Yn (boa te eibi, mV Aa na 
second group, 
He — —m40, +Å m; of, = — My Gn, + Aym,; 
PS 00, 
sv, A oV, ; 
NC —mjb tA,m; OU —m,b,, + Aam; 
SV, eV, , 
m. —m40€,4 T Agm,; MS — T, Cin + Aga; 
and finally, 
oV, 
ur 


(M*.) 


(N*.) 


(0*.) 


With respect to the six multipliers A,, Àz, Az, A,, A2, Ag, which were introduced by the 
3 final equations of condition (45.), and by the 3 analogous initial equations of condition, 


x.ma,=0, X.mb,=0, X.mc,=0; 
we have, by differentiating these conditions, 


x.mxzi=0, X.my;20, X.mz'=0, 


and 


2.ma,—-0, 2X.mb,-0, X.mc,=0; 


and therefore 


and 


eV, eV, èy, 
E rato 
En" ^ im' 1" Xm 
Eee um 
Be, gh Be 
XM a eee ANE 


(55.) 


(56.) 


(57.) 


(58.) 


(59.) 


11. As an example of the determination of these multipliers, we may suppose that the part 
V,, of the whole action V, has been expressed, before differentiation, as a function of H,, and 
of these other 6n — 6 independent quantities 


2,1—2,4,76), 

Va1—V,n7": 

2,17 9,n 77 bii 
and 

Q,1—0,,704, 

b —5,,—,, 

C,1—-C nm => 


Bn Tn Ca 1: 
Y, n-1 ~ Yin 7 1n-1» 
2 n-1  ^?,n ^ (MM 


AST PES OMS 
9,2 —9,n "l2; 
2,2 n = Cs, 


A, — 0,4 7X; 


b,2—-O n= Ba; 
€,2 — C,n — ya: 


Q,n—1 7 Un = X1; 
LINT T b = Bua , 
Cini on ^ Yn-1: 
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that is, of the differences only of the centrobaric coordinates*; or, in other words, as a function of 
the coordinates (initial and final) of » — 1 points of the system, referred to the ntt point, as an 
internal or moveable origin: because the centrobaric coordinates z,;, y,;, 2,;, %,;, 0,4, c,; may 
themselves, by the equations of condition, be expressed as functions of these, namely, 


x.m L.m L.m 
v,4,—6,— S. Y= mn Ser 2,4=%;- x. (62.) 


and in like manner, 


X.ma« L.m L.m 
a, = 04— =m’ b, B;— p C= Yi m (63.) 


in which we are to observe, that the six quantities én, Nns “ns «,, Bn» y, must be considered as 
separately vanishing. When V, has been thus expressed as a function of the centrobaric co- 
ordinates, involving their differences only, it will evidently satisfy the six partial differential 
equations,T 


yh, ^ 0, 5% <0, zo 

oa by, 8z, (Pi) 
V. òy, 8V, 
5a, ‘=0, aor at As. 


after this preparation, therefore, of the function V, , the six multipliers determined by (58.) and 
(59.) will vanish, so that we shall have 


^ = 0, Àz = 0, As = 0, A, = 0, Aa = 0, As = 0, (64.) 
and the groups (M!.) and (N!.) will reduce themselves to the two following: 
V, , eV, , ôy, ree 
Bm. a Bx. aum Sese ae i zv 
8V / dV, , eV, , 
Sy myn By, apa ose Sy =MnY ins (Q?.) 
, , "m 
òV, , + , eV, , 
Bron ELO e gg aon 
and 
9V, FS à L , èy, , 
õa —m40G,; Sa ai. — m4, Sa ^w —m OO, 
èy, p ee, ; oV, , 
8b i —mb; 5b. — mab;;; Mite vient 73 —m,b,,; (R+.) 
, ' "m 
òy, , 9V, fs èy, , 
$e. —m,0€,4; 8c E —WaQ€,; +++ 8c ^i —nT,n; 


analogous in all respects to the groups (C.) and (D.). We find, therefore, for the relative motion 
of a system about its own centre of gravity, equations of the same form as those which we had 
obtained before for the absolute motion of the same system of points in space. And we see that 
in investigating such relative motion only, it is useful to confine ourselves to the part V, of our 
* [The term centrobaric was first introduced by Móbius.] 
T E V, is expressed as a function of the é, n, £, «, B, y, 


SV, 8V, 
Stp 55 


($21, 2, ..., n= 1), 
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whole characteristic function, that is, to the relative action of the system, or accumulated living 
force of the motion about the centre of gravity; and to consider this part as the characteristic 
function of such relative motion, in a sense analogous to that which has been already explained. 


This relative action, or part V,, may, however, be otherwise expressed, and even in an 
infinite variety of ways, on account of the six equations of condition which connect the 6n 
centrobaric coordinates; and every different preparation of its form will give a different set of 
values for the six multipliers A,, Àa, A5, A,, Az, Ag. For example, we might eliminate, by a 
previous preparation, the six centrobaric coordinates of the point m, from the expression of V, , 
so as to make this expression involve only the centrobaric coordinates of the other » — 1 points 
of the system, and then we should have 


y SR AME NUNC a AND AM i 
ora i oce COM R. CM NG Gr) 


and therefore, by the six last equations of the groups (M1.) and (N1.), the multipliers would 
take the values 


À= — Tin, À= — Vn: A3= — Zin» A= n> Aga —b,,, Ag — 6,5; (65.) 


and would reduce, by (60.) and (61.), the preceding 6» — 6 equations of the same groups (M!.) 
and (N!.), to the forms 


V, BV, BV, 


is Et San, = Mabe» 2E manas 
gy. cmn y cms ym LT hihi (T1.) 
n Stemt, icc m m oom 
Ri — My % , mo =Ma, + oe e gins aah 
= -m fi, J= -Mfz ++ z= ii EN Bas (U?.) 
a —hy ki — Hays, ae —h yn 


12. We might also express the relative action V, , not as a function of the centrobaric, but 
of some other internal coordinates, or marks of relative position. We might, for instance, express 
it and its variation as functions of the 6n — 6 independent internal coordinates é, 7, £, «, B, y 
already mentioned, and of their variations, defining these without any reference to the centre 
of gravity, by the equations 

i= ti Ens M=VYi-Yn> MAN (66.) 
«,—0a,—0,, p,;—b;—b,, "y;—06;—6,. 
For all such transformations of dV, it is easy to establish a rule or law, which may be called the 
law of varying relative action (exactly analogous to the rule (B!.)), namely, the following: 


| ez 8, -z (552) Be, 2, o2. 84, 42-4, 30, (Vi.) 
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which implies that we are to express the half 7’, of the relative living force of the system as a 
function of the rates of increase 7’ of any marks of relative position; and after taking its variation 
with respect to these rates, to change their variations to the variations of the marks of position 
themselves; then to subtract the initial from the final value of the result, and to add the 
variations of the final and initial functions $, , 6, , which enter into the equations of condition, 
if any, of the form 4, — 0, ©, — 0, (connecting the final and initial marks of relative position,) 
multiplied respectively by undetermined factors A, , A,; and lastly, to equate the whole result 
to ôV, — 18H,, H, being the quantity independent of the time in the equation (50.) of relative - 
living force, and V, being the relative action, of which we desired to express the variation. It is 
not necessary to dwell here on the demonstration of this new rule (V!.), which may easily be 
deduced from the principles already laid down; or by the calculus of variations from the law of 
relative living force, combined with the differential equations of the second order of relative 
motion. 


But to give an example of its application, let us resume the problem already mentioned, 
namely to express ôV, by means of the 6n — 5 independent variations 93£;, ng, Òli, i, 58;, dy;, 
5H,. For this purpose we shall employ a known transformation of the relative living force 27’, , 
multiplied by the sum of the masses of the system, namely the following: 


2T,Xm-EX.mm,((z; — 24)? + (yi — Ye)? + (i — 2:5): (67.) 
the sign of summation X extending, in the second member, to all the combinations of points 
two by two, which can be formed without repetition. This transformation gives, by (66.), 

2T, Em-m,X,.m(£"?--3?--07?) 
+E, .mymy {l(i — E) + (ni — m)* + E-W) 
the sign of summation E, extending only to the first n— 1 points of the system. Applying, 


therefore, our general rule or law of varying relative action, and observing that the 6» — 6 
internal coordinates £, 7, C, «, B, y are independent, we find the following new expression: 


(68.) 


` 


SV, 2 18H, +=" . E, m (£8£ — a'a + 3/99 — BSB + L'8L —y'By) 


Xm 
gr E, mum, (£i E (BEBE) + (i ni) Gn) EG 00-969). | WE 
= sh. E, .mymy (ai — ai) (Bay — 803) + (Bi — Bi) Bi- Bx) + (71—74) y òy): 
which gives, besides the equation (O!.), the following groups: 
Seas Io4-0-m(u-E. 
E zi. E.m (ni =n) = m, (i- e. (X*.) 
BV, m, E mi 


BE. Em om =) =m (fi Xm 5 
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and 


eV, —m, 2 "Wn , à, m« 
7 Pai ^ Mia ) = -m (4 - žE), 
ôy, en — M; , "M , X, mp’ 1 
SLT mE -m (B-A), (Y*.) 
oV, —m; ; nah , X, my’ A 
greg dum oA manm ntis (mE) 
results which may be thus summed up: 
òV, — 18H, +E, . m (£'8£ — a'8a + $89 — B'8B + (86 — y'9y) 
1 ; , " 
— 5, 0m .,msE+X, mr’ . 3j, mòn - Xi, mU . Xi, mòt) (Zi) 
Fg, C, ma! E, mòa+ E,mf' .Z, mf +E, my E, my) 


and might have been otherwise deduced by our rule, from this other known transformation 
or T... 
(Z, më’)? + (E, mn’)? + (E, mb’)? 
22m ; 
And to obtain, with any set of internal or relative marks of position, the two partial differential 
equations which the characteristic function V, of relative motion must satisfy, and which offer 
(as we shall find) the chief means of discovering its form, namely, the equations analogous to 
those marked (F.) and (G.), we have only to eliminate the rates of increase of the marks of 
position of the system, which determine the final and initial components of the relative velocities 
of its points, by the law of varying relative action, from the final and initial expressions of the 
law of relative living force; namely, from the following equations: 


T,—$£E,.m(£? - 9? - 0?) — (69.) 


T,=U+H,, (50.) 
and 
T eX U,+4H,. (70.) 


The law of areas, or the property respecting rotation which was expressed by the partial 
differential equations (P.), will also always admit of being expressed in relative coordinates, 
and will assist in discovering the form of the characteristic function V,; by showing that this 
function involves only such internal coordinates (in number 6n — 9) as do not alter by any com- 
mon rotation of all points final and initial, round the centre of gravity, or round any other 
internal origin; that origin being treated as fixed, and the quantity H, as constant, in deter- 
mining the effects of this rotation. The general problem of dynamics, respecting the motions of 
a free system of n points attracting or repelling one another, is therefore reduced, in the last 
analysis, by the method of the present essay, to the research and differentiation of a function 
V,, depending on 6» — 9 internal or relative coordinates,* and on the quantity H,, and satis- 
fying a pair of partial differential equations of the first order and second degree; in integrating 
which equations, we are to observe, that at the assumed origin of the motion, namely at the 


* [The relative configuration of the 2(»— 1) initial and final positions of the planets is defined by 6» — 12 
mutual distances and in addition we have the three coordinates of any one of them relative to the Sun.] 
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moment when ¢=0, the final or variable coordinates are equal to their initial values, and the 


ô 

5H, 
differential alterations of the eoutdinates have ratios connected with the other partial differential 
coefficients of the characteristic function V,, by the law of varying relative action. It may be 
here observed, that, although the consideration of the point, called usually the centre of gravity, 
is very simply suggested by the process of the tenth number, yet this internal centre is even 
more simply indicated by our early corollaries from the law of varying action; which show that 
the components of relative final velocities, in any system of attracting or repelling points, may 


be expressed by the differences of quantities of the form bell Bhai i.d. : and therefore that 


m dx’ m dy’ m dz 
in calculating these relative velocities, it is advantageous to introduce the final sums Xm, Emy, 
Emz, and, for an analogous reason, the initial sums Ema, Emb, Ximc, among the marks of the 
extreme positions of the system, in the expression of the characteristic function V; because, in 
differentiating that expression for the calculation of relative velocities, those sums may be 
treated as constant. 


On Systems of two Points, in general ; Characteristic Function of the motion of any 
Binary System. 


13. To illustrate the foregoing principles, which extend to any free system of points, how- 
ever numerous, attracting or repelling one another, let us now consider, in particular, a system 
of two such points. For such a system, the known force-function U becomes, by (2.), 

U —m,msf (r), (71.) 
r being the mutual distance 
r— V (z, — 2a)? + (y, — Ya)? + (24 — Z2)? (72.) 
between the two points m, , Ma, and f (r) being a function of this distance such that its derivative 
or differential coefficient f’ (r) expresses the law of their repulsion or attraction, according as it 
is positive or negative. The known differential equations of motion, of the second order, are 
now, by (1.), comprised in the following formula: 


m, (x82, + yf dy, + 2; 024) + Mg (15917, + 1/5033 + 25025) = m, ma 9f (r); (73.) 
they are therefore, separately, 


f (r f (r : of (r 
Ms I ) ym TO i e; =M, To, 
of (r E r 2 òf (r 
fO. yam LO), gam » YO. 
The problem of integrating these equations consists in proposing to assign, by their means, 
six relations between the time ¢, the masses m4, Ma, the six varying coordinates %1, yi, 2,, 
X5, Ja, Z2, and their initial values and initial rates of increase a}, bi, c4, Qa, bg, Co, a1, 54, Cis 
a, b3, c5. If we knew these six final integrals, and combined them with the initial form of the 
law of living force, or of the known intermediate integral 


m, (x4? + yy? + 24) rm (x3? + Ya" +29") =m, mf (r) +H; (75.) 


r= 
(74.) 


£3 =M; 
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that is, with the following formula, 
im, (a5? + 03? + c3?) + 1m, (a5? + b;? + 02?) =m, ms f (r9) + H, (76.) 
in which 7, is the initial distance | 
ro — V (a, — a9)? + (b, — ba)? + (e — 63)’, (77.) 


and H is a constant quantity, introduced by integration; we could, by the combination of 
these seven relations, determine the time f, and the six initial components of velocity a}, bj, cj, 
a3, b3, c5, as functions of the twelve final and initial coordinates x,, Y1, 21, £a, y, 23, A, 54, €4, 
Mg, bg, c4, and of the quantity H, (involving also the masses:) we could therefore determine 
whatever else depends on the manner and time of motion of this system of two points, as a 
function of the same extreme coordinates and of the same quantity H. In particular, we could 
determine the action, or accumulated living force of the system, namely, 


t t 
Vom, | (e+ 2+ 218) dtm, [ (ag? typ? +252) a (A2) 
0 


as a function of those thirteen quantities x; , Y1, 21, 13, Yo, % 94, 01, Cy, Mp, by, Cg, H: and might 
then calculate the variation of this function, 


oV eV oV pa eV oV 
pest o 82, MEN St 924 Fir Sade 7s Syst Sz, O25 


— i o oV di oV eV 


e 
rr 5g 2 


86s (B?.) 


But the essence of our method consists in forming previously the expression of this variation, by 
our law of varying action, namely, 
8V =m, (2,92, — a; 8a, + y dy, — 019b, +21 92, — c480,) 
+ Mg (25 83, — 0580 + 150y5 — b59b, + 25824 — 05905) (C2.) 
+tdH; 
and in considering V as a characteristic function of the motion, from the form of which may be 


deduced all the intermediate and all the final integrals of the known differential equations, by 
resolving the expression (C?.) into the following separate groups, (included in (C.) and (D.),) 


oV oV t; OV 


Sa, = MX}, By, Mmh Se, 1 D) 
oV oV "n 
Ba, 7a By, mau? Sza = M22; 
and 

E e S MU E seas Inia 
8a, AU? 8b, 1*5: 8c, 1 e| "i 
AE 29 b! eK. ^ "i e 
Sa TEAN 8b, — M02, Ms —maes; 

HMPII 17 
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besides this other equation, which had occurred before, 
oV 


SH ‘ (E.) 


By this new method, the difficulty of integrating the six known equations of motion of the 
second order (74.) is reduced to the search and differentiation of a single function V; and to 
find the form of this function, we are to employ the following pair of partial differential equations 
of the first order: 


an (5) * i) "ms * Gs) ] mme. m 


ETET ET ab (EG QD] moon m 


combined with some simple considerations. And it easily results from the principles already laid 
down, that the integral of this pair of equations, adapted to the present question, is 


efe) m 


in which z, , y, , 2, , à, , b, , c, denote the coordinates, final and initial, of the centre of gravity 
of the system, 


V=V(c,=4,) + (y, —b,) + (e, e, .V 2H, (m, +m) + —— 


_ Matı + Mata ai a Ma TMs it! as Sh 

^ m 4- This of Mı + Mz ; m + ma (78.) 
_ 40, + Mz My _ mb, +mMgb, _ M,C, + OC, 
. tytn, ° “Onset” "> ee, 


and 3 is the angle between the final and initial distances r, ry: we have also put for abridgement 


p= + 2(m+ms(f ri J-5 


m,m,) r? (79.) 
the upper or the lower sign to be used, according as the distance r is increasing or decreasing; 
and have introduced three auxiliary quantities h, H,, H,, to be determined by this condition,* 


0-94 f zt dr, (I2.) 


combined with the two following, 


sat. gr dr Vs ca) rs Y Ge ae | (K2.) 


H,+H,=H; 


which auxiliary quantities, although in one view they are functions of the twelve extreme 
coordinates, are yet to be treated as constant in calculating the three definite integrals, or 
limits of sums of numerous small elements, 


òp p 
E pdr, L 55 d" T ur. dr. 


* [Cf. Appendix, Note 2, pp. 613-621.] 
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The form (H?.), for the characteristic function of a binary system, may be regarded as a 
central or radical relation, which includes the whole theory of the motion of such a system; so 
that all the details of this motion may be deduced from it by the application of our general 
method. But because the theory of binary systems has been brought to great perfection 
already, by the labours of former writers, it may suffice to give briefly here a few instances of 
such deduction. 


14. The form (H?.), for the characteristic function of a binary system, involves explicitly, 
when p is changed to its value (79.), the twelve quantities x, , y, , z, , 4, ,,, c, T, To, 9S, h, H,, H,, 
(besides the masses m,, m, which are always considered as given;) its variation may therefore 
be thus expressed: 

oV aV FLA eV oV oV 


Vds ea ys t &z, zc OF Ta da ot 3, òb, + Se, | òc, E 
58V SV sy 8V 8V ôy í 
tar +5, òro +55 89 tgz + 8H. SH +e 8H,. 


In this expression, if we put for abridgement 


" 2H, (m, +m) 
AP icis dm OU PY s, 08 e 


we shall have 


SK =A ty— dy), 3 =AYy—Dy 5 = ACEC) 
àz, òy, òz 
èy 8V V eg 
heehee age 8b, —A(b, —¥,); 8c, =A(c, -$,i 
and if we put 
Sye (reot -5 (81. 
Po= + ,[ 2(m, +m, ( (ro mm, 73’ . 
the sign of the radical being determined by the same rule as that of p, we shall have 
òV _ mmap eV — — m,fn»pg eV _ m; mah | (N?) 
ór m,4m, Sry mitm,’ 5% m +my ] 
besides, by the equations of condition (I?.), (K?.), we have 
oV 
ra 0?, 
s; (0%) 
and 
eV  8V  [*dr 
ai =$H. P*. 
. 8H, 78H.^ ot 6H,+6H,=5H i.) 


The expression (L?.) may therefore be thus transformed: 
ôV — A((z, —&,) (dz, "Ux 8a,) a (y, —b,) (8y, di b ») T (2, as s) (82, Ys 8e,)] 
rdr 


mum, — .oH; E 
tem (pòr — po9ro + hdd) + ay (Q?.) 
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and may be resolved by our general method into twelve separate expressions for the final and 
initial components of velocities, namely, 


A RERO UM Me or x] 
7m| mtm," o ata, Pie x," 
Lar 94 M ( Sr a 
{=——— =- " + — +h , 
m, y, m, tm, V } m, + ms Py By, 
gas ao clio quim (oi ^x. 
m, 52, M,+Ms My+ Ma \ Om — 82 (R2) 
aap da AM ai ~ )+ ir eri iz) 
a mass, mpm, " "| mum. V da, da)” 
1 ôy A m òr 99 
9—— = " b, + i ( hg). 
Mayo ™m,+m d m, +m, V By, Bye 
2 Papin w; [9r BS 
^1 ms 82; m+ m ") tmm, oi): 
and 
;» Mop allied, m, ( By d 
^17 m, 8a, mtm, + mem (ois "Sa, 
SE EEE ma f- ro a2) 
ow m rm, d e portare (ovs : h b) 
veo. ea. m, ( Bn 22) 
1 m, & m,-4 m n mm, Q c, B, ; (S?.) 
peter. SA (x jE eei ( ae | 
2 m, Sag My +My” My +My\Saq Bay)’ 
errem S at) 
uu mm, 2” lie (nit-^ : 
MERNE AV m, ( ro , 89). 
My 5o, m +m, * n) tmtm (rose x) 
besides the following expression for the time of motion of the system: 
8V "dr 
and 1 9) 4 a 2 
t=), s (T?) 
which gives by (K?.), and by (79.), (80.), 
tan Ru (U2.) 


The six equations (R?.) give the six intermediate integrals, and the six equations (S?.) give 
the six final integrals of the six known differential equations of motion (74.) for any binary 
system, if we eliminate or determine the three auxiliary quantities h, H,, H,, by the three 
conditions (I?.) (T?.) (U?.). Thus, if we observe that the distances r, ry, and the included angle 9, 
depend only on relative coordinates, which may be thus denoted, 

4,—2,4—$, Yr—-Yo=; xem 


(82.) 
d,—0a,—«*, b,—by—, c,—059—y, 
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we obtain by easy combinations the three following intermediate integrals for the centre of 
gravity of the system: 


x t=2,—Q,, y, t— y, —5,; z,t=2,—C,, (83.) 
and the three following final integrals, 
a t=, =a, bt=y,—b,, c,t=z,—C,, (84.) 


expressing the well-known law of the rectilinear and uniform motion of that centre. We obtain 
also the three following intermediate integrals for the relative motion of one point of the 
system about the other: 


EE d O 
* apt ar: 
AES. 2 BD 
4 BE. 0D 

and the three following final integrals, 
jo ro $09 
Pree cL 
R dry 99 
B —Pvgg hsp (86.) 
eni Org 99. 
Y-Ps ES 


in which the auxiliary quantities h, H, are to be determined by (I?.), (T?.), and in which the 
dependence ofr, 7), 3, on £, x, £, «, B, y, is expressed by the following equations: 


r=VE4 P, HI (87.) 
rro cos 9 — £a -- 9B + Cy. E 
If then we put, for abridgement, 
penpe VAT |, + ps voiMnoisergyo spi iod 
NDERIT YS - frySin 9^ TT Cina ied. 
we shall have these three intermediate integrals, 
éE = A£— Ba, n =An— BB, C’=Al— By, (89.) 
and these three final integrals, 
a’ = BE—Ca, B'—Bq—OB, y'=Bt- 0y, (90.) 
of the equations of relative motion. These integrals give, 
én! — në’ zu a ' Ti Box’ = B (an - Bé), 
nb’ — Uy = By’ —yB' = B (Bt — yn), (91.) 
LE — £0" — ya! — ay’ = B (y£— at), 
and 
C («B' -- Bau’) + £ (By’ —yB’) +» (ya! — ay’) 40; (92.) 
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they contain therefore the known law of equable description of areas, and the law of a plane 
relative orbit. If we take for simplicity this plane for the plane £x, the quantities £, £', y, y’ will 
vanish; and we may put, 


€=rcos0, | q—rsin6, P (93.) 
a — r9 COS Ó, B=rysin 4%, y=0, . 
and 
£-—r'cosÜ—ÜO'rsinÓ, ^ «'-—r'sin6-FÓ'rcos6, bra) (94.) 
a’ —rgcoS 0, — Oorosin fo,  B'—rosin0,--0,rocos0,,  y'—0, ' 


the angles 0, 0, being counted from some fixed line in the plane, and being such that their 
difference 


These values give 
£n, — né 2 1?0',  aB'— Ba’ =720,, an—BE=rrysin4, (96.) 
and therefore, by (88.) and (91.), 
1*ü' «120, — h; (97.) 


the quantity $4 is therefore the constant areal velocity in the relative motion of the system; a 
result which is easily seen to beindependent of the directions of the three rectangular coordinates. 
The same values, (93.), (94.), give 
£cos 0 -- sin 8 =r, £ cos 6+ 7'sin@=7", r pecia v Gs d (98.) 
«c080,4- Bsin0 —r,, «'cos0,-- B'sinOy—r,, £cos0,-- sin 0, —rcos9, ; 
and therefore, by the intermediate and final integrals, (89.), (90.), 


f'zp, fo po; (99.) 
results which evidently agree with the condition (T?.), and which give by (79.) and (81.), for all 
directions of coordinates, 

N. E fii 1). 
r T — 2 (m, +m) f(r)=ro ta ?Pmimijj(n)-? , m, m, ; (100.) 
the other auxiliary quantity H, is therefore also a constant, independent of the time, and enters 
as such into the constant part in the expression for ( 24 — = 3) the square of the relative velocity. 


The equation of condition (I?.), connecting these two constants h, H, , with the extreme lengths 
of the radius vector r, and with the angle 9 described by this radius in revolving from its initial 
to its final direction, is the equation of the plane relative orbit; and the other equation of 
condition (T?.), connecting the same two constants with the same extreme distances and with 
the time, gives the law of the velocity of mutual approach or recess. 


We may remark that the part V, of the whole characteristic function V, which represents 
the relative action and determines the relative motion in the system, namely, 


=f Ee 1 2 
y - à x (io +f par), (V*.) 
may be put, by (I?.), under the form 
mima [* àp i 
V s : |J. (»- ^ai) dr, (W?.) 
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or finally, by (79.), 


V,=2 tS Oe ede; (X2.) 
the condition (I?.) may also itself be transformed, by (79.), as follows: 
wol v iie A (Y2.) 
To rp p 


results which all admit of easy verifications. The partial differential equations connected with 
the law of relative living force, which the characteristic function V, of relative motion must 
satisfy, may be put under the following forms: 


eV, |? dV, 2m, Mz 
(=) tals) = m, +m ha | 


dV, \2 òy, 2m,m; 
EARETCO BE (Cot 


and if the first of the equations of this pair have its variation taken with respect to r and 9, 
attention being paid to the dynamical meanings of the coefficients of the characteristic function, 
it will conduct (as in former instances) to the known differential equations of motion of the 
second order. 


(Z?.) 


On the undisturbed Motion of a Planet or Comet about the Sun : Dependence of the 
Characteristic Function of such Motion, on the chord and the sum of the Radii. 


15. To particularize still further, let 
1 - 
fr) r? (101.) 


that is, let us consider a binary system, such as a planet or comet and the sun, with the New- 
tonian law of attraction; and let us put, for abridgement, 


2 


Mı T TR — p, "is 2H, =a. (102.) 
The characteristic function V, of relative motion may now be expressed as follows: 
Iu T 3 TTIE. ) A3 
ya (ovo [+ J2-1-5 ar); (A) 


in which p is to be considered as a. function of the extreme radii vectores r, ry, and of their 
included angle 9, involving also the quantity a, or the connected quantity H, , and determined 
by the condition 


9 — ———————— , B?. 
j d: P AC yc 
» — — — — — 
" rp ap r 
that is, by the derivative of the formula (A?.), taken with respect to p: the upper sign being 
taken in each expression when the distance r is increasing, and the lower sign when that distance 
is diminishing, and the quantity p being treated as constant in calculating the two definite 


integrals. It results from the foregoing remarks, that this quantity p is constant also in the 
sense of being independent of the time, so as not to vary in the course of the motion; and that 
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the condition (B?.), connecting this constant with r, 79, 9, a, is the equation of the plane relative 
orbit; which is therefore (as it has long been known to be) an ellipse, hyperbola, or parabola, 
according as the constant a is positive, negative, or zero, the origin of r being always a focus of 
the curve, and p being the semiparameter. It results also, en the time of motion may be thus 
ex pressed: 


t=—+ =——_ , (C3.) 
and therefore thus: | 
t= 3 (D?.) 


which latter is a known expression. Confining ourselves at present to the case a> 0, and intro- 
ducing the known auxiliary quantities called excentricity and excentric anomaly, namely, 


i fict (103.) 


v=cos-} end b j (104.) 


and 


which give 
t V/2ar — 1? — pa — aesin v, (105.) 
v being considered as continually increasing with the time; and therefore, as is well known, 


r=a(l—ecosv), ry=a(l—ecosyp), 


1 106. 
9 —2tan^ x “tanz — 2tan^1 h SN tan 2 j ( ) 
1—e 2 
and 
a? D J 
t= /—.(v—vp—esinu+esin vo); (107.) 
p 


we find that this expression for the characteristic function of relative motion, 


r /2 1 
i -i)a 
V aae ma p AE aait aA eae) 


deduced from (A?.) and (B*.), may be transformed as follows: 


Y, «mm, J% (os eine —esinuy); f (F?.) 
in which the excentricity e, and the final and initial excentric anomalies v, vg, are to be con- 
sidered as functions of the final and initial radii r, ry, and of the included angle 9, determined 
by the equations (106.). The expression (F?.) may be thus written: 


V, e NE (v, 4- e, sin v,), (G?.) 
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if we put, for abridgement, 

ui Prog ggg OE. 
, 2 ? , 9 > 
for the complete determination of the characteristic function of the present relative motion, it 
remains therefore to determine the two variables v, and e, , as functions of r, ry, 9, orof some other 
set of quantities which mark the shape and size of the plane triangle bounded by the final and 
initial elliptic radii vectores and by the elliptic chord. 


(108.) 


For this purpose it is convenient to introduce this elliptic chord itself, which we shall call 
+7, so that 
T? =r? + r? — 2rro cos Ñ; (109.) 
because this chord may be expressed as a function of the two variables v, , e,, (involving also 
the mean distance a,) as follows. The value (106.) for the angle 9, that is, by (95.), for 0 — 0,, gives 


g—2tan-|, JFS tan; = & - 2tan-i| 1+ tan zj- o, (110.) 


w being a new constant boi E of the time, namely, one of the values of the polar angle 6, 
which correspond to the minimum of radius vector; and therefore, by (106.), 


rcos(9—a@)=a(cosv—e), rsin(0—c)—a'V1-—e?sinv, | (111) 
ro eos (0 — w) — a (cosvọ— e), rosin (0, — w) 2 a'V1— e?sin v; 
expressions which give the following value for the square of the elliptic chord: 
= {r cos (0 — w) — ra cos (05 — w)} + {r sin (0 — w) — ry sin (0) — &)) 
= a? {(cos v — cos vo)? + (1 — e?) (sin v — sin v,)?) 
= 4a? sin v? {(sin Z ye) +(1—e?) (cos id 29 yr 


= 4a? (1 — e?) sin v?: 


we may also consider 7 as having the same sign with sinv,, if we consider it as alternately 
positive and negative, in the successive elliptic periods or revolutions, beginning with the 
initial position. 

Besides, if we denote by o the sum of the two elliptic radii vectores, final and initial, so that 


g —T Te, (113.) 
we shall have, with our present abridgements, 
o = 2a (1— e, cos v,); (114.) 


the variables v, , e, are Pts "a. functions of c, 7, a, and consequently the characteristic function 
V, is itself à function of those three quantities. We may therefore put 
| Pa Tmaw 

m4 4- m, 

w being a function of o, 7, a, „of which the form is to be determined by eliminating v, , e, between 
the three equations, 


(H*.) 


pam 


w — 2'V ua (v, 4- e, sin v,), 
o= 2a (1— e, cos v,), (D.) 
7— 2a (1— e?) sin v,; 


HMPII 18 


www.rcin.org.pl 


138 II. GENERAL METHOD IN DYNAMICS ) [15 


and we may consider this new function w as itself a characteristic function of elliptic motion; 
the law of its variation being expressed as follows, in the notation of the present essay*: 
dw = ESE — ada -- 5/93 — B'8B + C'96 — y'Sy + ien : (K?.) 
In this expression, é, 7, C are the relative coordinates of the point m, , at the time t, referred to the 
other attracting point m, as an origin, and to any three rectangular axes; £', n’, £' are their rates 
of increase, or the three rectangular components of final relative velocity; «, B, y, «', B’, y’ are the 
initial values, or values at the time zero, of these relative coordinates and components of relative 
velocity; a is a quantity independent of the time, namely, the mean distance of the two points 
m4, Ma; and p is the sum of their masses. And all the properties of the undisturbed elliptic 
motion of a planet or comet about the sun may be deduced in a new way, from the simplified 
characteristic function w, by comparing its variation (K?.) with the following other form, 
dw dw dw 


ass pit pte 8 3 
w= ôo + S ral da; (L2.) 


in which we are to observe that 
o= VE ry OY Va? + B+ y?, | (M3.) 
T= £V (£—a)** (o— B * (L—-yP*. | 
By this comparison we are brought back to the general integral equations of the relative 
motion of a binary system, (89.) and (90.); but we have now the following particular values for 
the coefficients A, B, C: 


18w 198w 18w 1w 18w. 4 
= bo. We DNE SEM Lia a oes e 
dw dw dw 


and with respect to the three partial differential coefficients $5 Ba? we have the following 


? ôo , 
relation between them: 

Ped Seri a (0?.) 

Sa 0o 87 2' i 
the function w being homogeneous of the dimension 3 with respect to the three quantities 
a, c, 7; we have also, by (I?.), 
e 1t sin v, i v1-e (P3.) 
a 


8o e,—cogsv, òr Na cor'v,—e,' 
ww  —2yur Ow (dw\? pw 4uo l n 
E CAT (Q3 


a or? 
from which may be deduced the following remarkable expressions: 


and therefore 


00 òr ott & (R2) 
(z-z)- dt. odi 
S80 87) c—-T a 


These expressions will be found to be important in the application of the present method to 
the theory of elliptic motion. 


* | aH, =a | 
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16. We shall not enter, on this occasion, into any details of such application; but we may 
remark, that the circumstance of the characteristic function involving only the elliptic chord 
and the sum of the extreme radii, (besides the mean distance and the sum of the masses,) 
affords, by our general method, a new proof of the well-known theorem that the elliptic time 
also depends on the same chord and sum of radii; and gives a new expression for the law of this 
dependence, namely,* 


esi, (S3.) 


We may remark also, that the same form of the characteristic function of elliptic motion con- 
ducts, by our general method, to the following curious, but not novel property, of the ellipse, 
that if any two tangents be drawn to such a curve, from any common point outside, these 
tangents subtend equal angles at one focus; they subtend also equal angles at the other. 
Reciprocally, if any plane curve possess this property, when referred to a fixed point in its own 
plane, which may be taken as the origin of polar coordinates r, 0, the curve must satisfy the 
following equation in mixed differences: 


A0 1 dl 
ootan (F) Az - (8-2) 557. (115.) 
which may be brought to the following form, : 
d | d*|l 
(p+); -* (116.) 
and therefore gives, by integration, 
P (117.) 


dex --ecos (0 — v): 

the curve is, consequently, a conic section, and the fixed point is one of its foci. 
The properties of parabolic are included as limiting cases in those of elliptic motion, and 

may be deduced from them by making 

H,=0, or &-oo; (118.) 
and therefore the characteristic function w and the time f, in parabolic as well as in elliptic 
motion, are functions of the chord and of the sum of the radii. By thus making a infinite in the 
foregoing expressions, we find, for parabolic motion, the partial differential equations 


dw Ow]? 4u w OwY 4u, : 
(ats TIFT Gana) V git b 
and in fact the parabolic form of the simplified characteristic function w may easily be shown 
to be 
w=2V u (Va 4-7 X Va —*), (U3.) 


T being, as before, the chord, and c the sum of the radii; while the analogous limit of the ex- 
pression (S?.), for the time, is 
tiec ua oe Sh (V3.) 
6v H 
which latter is a known expression. t 
* [Equation (S?.) is easily reducible to the well-known form of Lambert’s theorem. Cf. Whittaker, Analytical 


Dynamics (1927), p. 92.] 
T [This special case of Lambert’s theorem was given by Euler in 1742. Euler, Miscell. Berolin. Tome vr.] 


18-2 
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The formule (K?.) and (L?.), to the comparison of which we have reduced the study of 
elliptic motion, extend to hyperbolic motion also; and in any binary system, with Newton's 
law of attraction, the simplified characteristic function w may be expressed by the definite 


integral* 
ae Ws 
»-[ EE 4a T^ sd 


this function w being still connected with the relative action V, by the equation (H?.); while 
the time t, which may always be deduced from this function, by the law of varying action, is 
represented by this other connected integral, 


alf odum s ET $ 
=i (& £) n (X3.) 


provided that, within the extent of these integrations, the radical does not vanish nor become 
infinite. When this condition is not satisfied, we may still express the simplified characteristic 
function w, and the time ¢, by the following analogous integrals: 


w= |" = Leg ke | (Y?.) 
and 
ie [^s PM | (Z?.) 
in which we have put for abridgement 
o, 5, v eT, (119.) 


and in which it is easy to determine the signs of the radicals. But to treat fully of these various 
transformations would carry us too far at present, for it is time to consider the properties of 
systems with more points than two. , 


On Systems of three Points, in general; and on their Characteristic Functions. 


17. For any system of three points, the known differential equations of motion of the 2nd 
order are included in the following formula: 


m; (2,82, + f Sy + 25923) + ma (2592, + Y3 Ya +24 529) | — (120.) 
+ Mg (235925 + 5 9s +23 O23) = SU, 


the known force-function U having the form 
U =m m f2 +m, mf) + mmf, (121.) 
in which f(9, f(,9. f(,3 are functions respectively of the three following mutual distances of 
the points of the system: 
10,9 =V (y, —23)* + (y1 — Yo)" + (44 — 2a)", 
19 VG Gr t (s a (122) 
19,3 =V (arg — T3)? + (ya — Y3)? + (2a — 23)”: 


* [Cf. p. 48.] 
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the known differential equations of motion are therefore, separately, for the point m,, 
80.9 8fa.9 
"s. x US As y | 
Sfo 2 sf) 
Y= ms^.—— ghee sy (123.) 
»_ Sf 8f.» 
nae Ore tg eg ih 


with six other analogous equations for the points m, and ma; xi, &c., denoting the component 
accelerations of the three points m,, Ma, Mg, or the second differential coefficients of their co- 
ordinates, taken with respect to the time. To integrate these equations is to assign, by their 
means, nine relations between the time t, the three masses m,, Ma, Mg, the nine varying co- 
ordinates 21, Y1, 24, V2, Ya, 23; V3, Y3, 23, and their nine initial values and nine initial rates of 
increase, which may be thus denoted, a,, b1, C1, 04, by, Ca, Ag, bg, C3, à, Ois Cis Ag, 05, Cos Ag, 
bz, cz. The known intermediate integral containing the law of living force, namely, 


jm, (c3? + y1?* + 2;*) + $m (23? + YQ? + 2;?) + $m, (23? + ys? + 25") (124.) 
—mqimsf9?»--m,msf(*9--m,m;f9-- A, ; 
gives the following initial relation: 
$m, (a4? + by? +03”) + mg (a5? + b;? + 03?) + 3m; (a5? + 05? + c3?) (125.) 
. =m, Myf fh? +m m fi? + mm f? +H, ; 


in which fí^?, f(^9, f9 are composed of the initial coordinates, in the same manner as 
f (9, f(.9, f> are composed of the final coordinates. If then we knew the nine final integrals of 
the equations of motion of this ternary system, and combined them with the initial form (125.) 
of the law of living force, we should have ten relations to determine the ten quantities t, a; , b; , c1, 
3, 5, 65,05, 05, c5, namely, the time and the nine initial components of the velocities of the three 
points, as functions of the nine final and nine initial coordinates, and of the quantity H, in- 
volving also the masses; we could therefore determine whatever else depends on the manner 
and time of motion of the system, from its initial to its final position, as a function of the same 
extreme coordinates, and of H. In particular, we could determine the action V, or the accumu- 
lated living force of the system, namely, 


t t t 
Vom, | (+y ditm] (etyde | (er+ust+e)de (AS) 


as a function of these nineteen quantities, %1, Y1, 74; La, Yo, Zas Tg» Yz» 23» A, 01; Cy, Aq, Dg, Ce, 
45, 05, c4, H; and might then calculate the variation of this function, 


ax xo RMLCURE MERCI NINL A 
BV = ag, Mmi Mae bat ga, bat gp, beso e, 


nii aV SV àv as 8V 


ima e T 8yg + Tar hs Bb, = 90g at 5, 865 - 
.9V oV oV SV sa Na BV 
oV 
tag 9H 
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But the law of varying action gives, previously, the following expression for this variation: 
òV = m, (2182, — a; 8a, + y19y4 — 019b, + 2, 02, — c190,) 
+ Mg (2592, — 2590 + 59) — bz ba + 2592, — 05903) 
+ Mg (2582 — 25905 + 59 — bg 8b, + 25925 — 03903) 

+tdH; 
and shows, therefore, that the research of all the intermediate and all the final integral equations, 
of motion of the system, may be reduced, reciprocally, to the search and differentiation of this 


one characteristic function V; because if we knew this one function, we should have the nine 
intermediate integrals of the known differential equations, under the forms 


oV PER A 


(C*.) 


Sa, acis ày, tt oo 

oV " aV v , eV — , 4 
Bm, 0??? Sy, 242 82,0 2*2? (D*.) 
oV (OOK. x AV j 


Hai i Sy, 348 Bc ve 


and the nine final integrals under the forms 


NV oer CBR is Mew samen Do qu 
Sa, 1^1 8b, jf ©) Se) | due Ea 
8V d^ lia bed à; " 
Ba, Mata» Bb, aba» il^ — M204, (E4.) 
BY icosts. h as ae pae b; Ys — fli, C; 
Bag oia sto dili ca, d ol Asaa 


the auxiliary constant H being to be eliminated, and the time ¢ introduced, by this other 
equation, which has often occurred in this essay, 


t=—. (E.) 


The same law of varying action suggests also a method of investigating the form of this 
characteristic function V, not requiring the previous integration of the known equations of 
motion; namely, the integration of a pair of partial differential equations connected with the 
law of living force; which are 


Sn (ie) * (5g) * (5a) |g ae) * C7 eT | ga 
d PANA: EE a) | = m maf 0,9 + mms f+ m,m3f 99 + H, 
9] «8 G+) bei 


+f 


EAG] 


(5. 5) + *( =) = m mfi m,máff9 + mgm 6°) +H. 
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And to diminish the difficulty of thus determining the function V, which depends on 18 co- 
ordinates, we may separate it, by principles already explained, into a part V, depending only 
on the motion of the centre of gravity of the system, and determined by the formula (H!1.), 
and another part V, , depending only on the relative motions of the points of the system about 
this internal centre, and equal to the accumulated living force, connected with this relative 
motion only. In this manner the difficulty is reduced to determining the relative action V,; 
and if we introduce the relative coordinates 


£—X,—-X. 71—91— Vs; apd Sg (126.) 
És =£a— T3, Ne=Yo—-Ys, [37 73— 23, ; 

and %=A,—a3, P,=b,—bs, -— (127.) 
&370,—405, a—b3—05, ya—03— 05, 


we easily find, by the principles of the tenth and following numbers, that the function V, may 
be considered as depending only on these relative coordinates, and on a quantity H, analogous 
to H (besides the masses of the system); and that it must satisfy two partial differential equa- 
tions, analogous to (F*.) and (G*.), namely, 
(s; 
9b, 


im, se) 55). * (5) | ms se) 6) + i | 
1 (z xy (E y d Ay (H4) 


ocdsztszz] H sy 
2m, (A55 86) Vom 85) ^ GS Bd 


=m, mf? +m, msf09--m,msf*9?-- H,; | 
am; s) * (spi) + (5) | Fs) * 5.) 6) | 
2m, (Vo, 9B, (5; 2m (Aot, ps Sya 


dl òV, oV, "(se 52) "(t| (I*.) 
a, (oe i) 9B, Op, by, Sy, 


=m, Mmf? + my Ms f? +m m fH: 


and 


the law of the variation of this function being, by (Z?.), 

òV, —t8H , +m, (£19£, — 0190 +181 — By SB, + 0905, — y1934)] 
+ My (£588, — Aa Òa + 15903 — B39 Bs + 0,00, — yo 572) | 

(m4 £i + My £;) (m4 SE, + m48£,) — (m, o + ma a5) (My 9o, + Mo Òda) 
+ (m4 + Mana) (M4 ÒN + MaN) — (m4 By + ms B3) (m, 8B, + m38;) | 
+ (m4 C + ma 55) (m4 86, + m480,) — (m, y3 + may) (m, Sy, + Maya) J 
which resolves itself in the same manner as before into the six intermediate and six final 

integrals of relative motion, namely, into the following equations: 
1 V, ,  mi&-ms&. uL ux Mm & + MoE. . 
m, dé, ^^ my+mMy+mMs Md, 2 m +m, +m, 
1 àV —, mQqtm.m, 19V, ,_ mni tman, 
my 5, 1 my +m+my MySy "My +M_y+M, 
DN, OQ mitt, iN, Q mt tml 
mob, ! m.m. mba 75 m,+m,+mg 


(K*.) 


Mı + Ma + Mg 


(L*.) 
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and «e UV, , My oc + Mg My —18V, =o! _ matma, 


m, Sa, mQcmO-Rma My St, My +My +Ms’ 
2185 Sp. ant NODE y. AE ets iU UM. (M4) 
m, 9p, ki M+M, +M; ms SB, £ m. + Ta + Mg’ 


-18V, __,_myi+Mmays, —18V, __, _ myi +My, 


m, y, - mı +M +m My dy, ** mQ4mag4ms 


DL (01.) 


18. The quantity H, in V,, and the analogous quantity H, in V, , are indeed independent 
of the time, and do not vary in the course of the motion; but it is required by the spirit of our 
method, that in deducing the absolute action or original characteristic function V from the 
two parts V, and V, , we should consider these two parts H, and H, of the original quantity H, 
as functions involving each the nine initial and nine final coordinates of the points of the 
ternary system; the forms of these two functions, of the eighteen coordinates and of H, being 
determined by the two conditions, 


ov, Or, re 4 
BH, 8H. H,+H,=H. (N*.) 


However, it results from these conditions, that in taking the variation of the whole original 
function V, of the first order, with respect to the eighteen coordinates, we may treat the two 
auxiliary quantities H , and H, as constant; and therefore that we have the following expressions 
for the partial differential coefficients of the first order of V, taken with respect to the coordinates 
parallel to z, 


MO, m MN, NS. 0m M, 

Sa, Sé, m,+m,+m, ÒL,’ da, Sa, m,+m,+mz, da,’ 

a A LAT 1 alte. sun OY PN 
Sé, mi +m, + mgr,’ Sa, Sa, Mı +M -+ M 54,” E 

» V, By m S, sov  &, ®, m, 003, 


Sr, BE, Òf, m,4m,vm,0m, Say — 0m, Say m, m.4-m;5a,' 


together with analogous expressions for the partial differential coefficients of the same order, 
taken with respect to the other coordinates. Substituting these expressions in the equations of 
the form (O.), namely, in the following, 


oV 8V 85V 38V BV 95V 


Bi, TRUE Sa TRUE 

oV èV 98 Aad eV əy 

uu L7 5b, ^ 85, b, 8b, =0, pt. 
By, By, By, 8b, vili 
eV èV èy ni o mode a 


id Ghee he, | Ate Ate 
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we find that these equations become identical, because 

eV, V, eV, V, 
Sx, Ba, y, 8b, - 

But substituting, in like manner, the expressions (O*.) in the equations of the form (P.), of 


which the first is, for a ternary system, 


E Ip NC UN NIC 
1 Sy, 9h ór,  ?*89, Y2 Sm, 7s ys Js 823 


aV, 


òV, 
x due 


mesh 8c 


0, t 


(Q*.) 


" " 


4 
a a ad a SLAP LP zm m i 
"A OO HET TA. (UN T 
and observing that we have 
oY, Vea or, —b S, 0, (84) 


E PA A A 

along with two other analogous conditions, we find that the part V,, or the characteristic 
function of relative motion of the ternary system, must satisfy the three following conditions, 
involving its partial differential coefficients of the first order and in the first degree, 


62 PRAM OM, EN A EA A 
1 Sm, "hag. 2 Sna UEA 1 9f, 18a, 2 5B. 2 8a," 
èy ôV, SV, òV èV èV èy, V, 
Ome the i ee a i du URB eA T4. 
7h 85 5 Sn, Ne JA b» Sno By 9y, i Y1 5B, B» y, Y2 9B, ( ) 


urs AI LAENE AN LADO leaky Heb gle 
X He) e M Obs %1 Sa, "y Lm 72 Sy. 


which show that this function can depend only on the shape and size of a pentagon, not generally 
plane, formed by the point m, considered as fixed, and by the initial and final positions of the 
other two points m, and ms; for example, the pentagon, of which the corners are, in order, 
Ms; (m4), (Ma), m, , Mı; (m4) and (ma) denoting the initial positions of the points m, and m, , referred 
to mz as a fixed origin. The shape and size of this pentagon may be determined by the ten 
mutual distances of its five points, that is, by the five sides and five diagonals, which may be 
thus denoted: 


M3 (m)=V58;, (m) (m3) — V 8s, (ma) ma — V ss, mam. — V 8,, mima =V 8s, 


128. 
Ms (m4) - V dy, (m4) mg — V d;, (m3) m, — V d, mam — V ds, Mı aaea À 


the values of sı, ... d; as functions of the twelve relative coordinates being 


sı=41 + Bi yl, Sa—(x3—04)*- (B — B)? + (ya — y)"; 
83 = (£a — 3)? + (qa — P2)? + (t2 — Y2)°, 

85— 64 mb +i, 84— (£1 — E)? + (m — 99)? + (5 — 53), (129.) 
d,=a5+ Bit yi, dy = (£4 — %1)? + (my — Bx)? t (E — y), 
d; = (£1 — ag)? + (m — Ba)? + (51 — v2)”, 

d,—£&it- «91-0, ds=(E,—%)? + (m — B)? + (5 79. 


HMPII I9 
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These ten distances Vs,, &c., are not, however, all independent, but are connected by one 


equation of condition, namely,* 


025151 +8383 + 82.8? 
+d? o 44 — 4d 
-didi — 4didà — edid 
— 283838}  —28528,84  — 282854, 
—25825,d,  —2s52s,d, —2s2s,d, 
—2s?8,d,  —2s528,d, | —25$5,d; 
—2s5,d,di —2s,d;d} —2s8,d,d? 
— 2s dåd,  —25,d2d;  —254did, 
— 2d dåd} —2d,d;d, —2d;d}d; 


+838? + 82.83 
+sidi +s3d3 
+ did; +d3d} 
—28?8,8, | — 2828,85 
—2s$s,d, | —2825,d, 
—252s,d, —2s?s,d, 


—28d,dj  —25,d,d; 

— 28, ddy —2s,d3d, 
—2d,djd, —2d;did, 
—48,8598,04 —485,8,85d4, —485858,d; — 48484890, — 4858585. 
—48,d, dd, —48,d,d,d; —48,d,d,d, —48,d,d,d, —48,d,dgd, 
—28,8,8,d, —28,838,d; —28,8,8,d, — 28485812 — 2858, 8,d3 

— 28,84d,d, —28,8,d,d, —28,8,d,d, — 2848,d,d, —28,8,d;d, 

— 2s d dd; —28,d_44d, — 255d5d,d, — 25,d, d, d, — 25d, d, d, (130.) 
+ 28182858, +282838485 +28384858, +284858)82 +2858, 8283 
+2818 84dg +2828385d4 --2555,8, dg + 28485890, +2858, 83d, 
4-28,858,d, +2828485da +283858,d, -- 28,8185 d, + 28582830; 
4-25,85d4d, + 2855,d,d; +28,8,d,d, + 28,85d, d, + 285,8, d, d. 

+ 25,85d, d, +28,8,dgd, -- 28455d,d; -- 25,8, d5d, + 28585 d, d. 
4-25,5,d, d, - 28585d, d, -- 2555, d, d, +2848 .dqd, + 2858505 d, 
4-25,5,d, d 4- 284585d, d, -- 2555, d5d; -- 25,85 d, d, + 28585d5 d. 
4-25,5,d,5d, -- 28585d4 d, -- 2855, d,d; + 28,8 .d5d, + 28585 d, d, 
-25,8,d5 d, -- 28585 d, d, -- 2848, dd, -- 28,85d, d, + 28555d, d. 
4-25, d, d, d, + 25,d, d, d, + 25d, d, d; + 28,d,d5 d, + 28545 d, d 

+ 25, ddd; + 25,d, d, d, -- 254d5d, d, + 28,d, do d, + 285d, dd, 

+ 2d, d, d, d, + 2d, d, d, d, + 2d, d, d d, + 2d, d; d, d, + 2d5d, d, d; J 


they may therefore be expressed as functions of nine independent quantities; for example, of 
four lines and five angles, (9, (9, r2, r(9, 0, 0P, 02, 019, ı, on which they depend as follows: 


* [This is the relation connecting the distances of five points in space. It can be written more symmetrically 
in determinant form (cf. Salmon, Geometry of three Dimensions (1914), p. 47). 


- 


* 


~ 


0 
1 
1 
1 
1 
1 


i L 

» Sa, 
&, 9, 
d,, ds, 
da, 8$, 
&, d, 
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8,7 79?, 

85 — r0? + r (22 — 2r Pr? (cos 0P cos 019 + sin 0j? sin 6$? cos), 
84 = 102 p (2 — 2/905 (9 cos (000 — 92), 

84 72 4 4002. 272) 2 (cos (X? cos 09 + sin (? sin 09 cos c), 
8,1002, 

d, —r$2?, 

dy — 102 4 y (02 — 2:9 (D (cos 99? cos 0(D + sin 0? sin 0j? cost), 
d, — r(92 (002 — X (7M (cos 0 cos 0V + sin 0f? sin 6 cos), 
d, — 192, 

d; — 1 (02 + p (02 — Vr (D cos (0V — 8(D), 


(131). 


the two line-symbols 7, 7 denoting, for abridgement, the same two final radii vectores which 
were before denoted by 759, r9, and r(9, r(? representing the initial values of these radii; 
while 00, 0@, 0(D, 0(? are angles made by these four radii, with the line of intersection of the two 
planes r(? r®, (9 59: and vis the inclination of those two planes to each other. We may therefore 
consider the characteristic function V, of relative motion, for any ternary system, as depending 
only on these latter lines and angles, along with the quantity H,. 


The reasoning which it has been thought useful to develope here, for any system of three 
points, attracting or repelling one another according to any functions of their distances, was 
alluded to, under a more general form, in the twelfth number of this essay; and shows, for 
example, that the characteristic function of relative motion in a system of four such points, 
depends on the shape and size of a heptagon, and therefore only on the mutual distances of its 


E : 7x6 
seven corners, which are in number D 


= | 21, but are connected by six equations of condition, 


leaving only fifteen independent. It is easy to extend these remarks to any multiple system. 


General method of improving an approximate expression for the Characteristic 
Function of motion of a System in any Problem of Dynamics. 


19. The partial differential equation (F.), which the characteristic function V must satisfy, 
in every dynamical question, may receive some useful general transformations, by the separa- 
tion of this function V into any two parts 


V,+V,=V. (U*.) 
For if we establish, for abridgement, the two following equations of definition, 
2 2 2 
Y MM 
2m (Vox by òz (V4.) 
Oe etek. a T (oV | 
Ti=2-am e) * Gv) * e] | 
analogous to the relation 
1 ((6VyY (8V|*? [(8V]? 
a£ cT via a Wt. 
7-2- ss) av) * e) | Be 
19-2 
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which served to transform the law of living force into the partial differential equation (F.); we 
shall have, by (U*.), 
òV, eV, dV, 8V, 98V, dV, 
TzT 
iT alse Sar t By By t Bz 3 
and this expression may be further transformed by the help of the formula (C.), or by the law 
of varying action. For that law gives the following symbolic equation, 
NANI Tt 
m\ da da" dy Sy 2902] dt 
the symbols in both members being prefixed to any one function of the varying coordinates of 
a system, not expressly aeo the time; it gives therefore by (U*.), (V*.), 
òV, òV, 3V,8V, 8V,9V;| dV, á 
Hr ôx by ws dz 1) dt Bie ut 
In this manner we find the following general and rigorous transformation of the equation (F.), 
dV, 
qi T- Toe Ty (A5.) 
T being here retained for the sake of symmetry and conciseness, instead of the equal expression 
U +H. And if we suppose, as we may, that the part V, , like the whole function V, is chosen so 
as to vanish with the time, then the other part V, will also have that property, and may be 
expressed by the definite integral, 


E j r- T, + T.) dt. (BS.) 


(X*.) 


(Y*.) 


More generally, if we employ the principles of the seventh number, and introduce any 3» 
marks 7 , 72, ... Ngn, of the varying positions of the n points of any system, (whether they be the 
rectangular coordinates themselves, or any functions of them,) we shall have 


eV èy SV 
T= rS bd A j C». 
omy” S12 Snan WE 
and may establish by analogy the two following equations of definition, 
r-r (1,3, 5), 
85, 99g! — 97g, (D5.) 
T,-F pa Vv, — oW | f 
° 82’ 7 Stan)’ 


the function F being always rational and integer, and homogeneous of the second dimension; 


and being therefore such that (besides other properties) 
9T, 9V, 57, SV, 9T, SV, 


Pa T T3 MI ect—RRi— ES. 
at Pat 57, 8m EA 5 A Bm at 
9m "M 99s 
07 | Od me oT, oT 9T, 135 (FS.) 
mds NA 5 ole 5 OV 5 YA ; 
8n 09 ôn 9935 nsn San 
and 
STs 8V, 8T, 9V, 9T, dV, 
—L 3 2 Bs’ —2T,. G5, 
~8V, 5 39, 8. "1 Ve Bee, a M 
071 na - ns» 


www.rcin.org.pl 


19] II. GENERAL METHOD IN DYNAMICS 149 


By the principles of the eighth number, we have also, 
Mv ee wp oT 


— — ——— — pete EMEN T ES 5 
R14 Ni» J4 a> «+s Ha N3n> (H°.) 
ôn 875 873, 
and since the meanings of 7, ... Ngn give evidently the symbolical equation, 
ò 6a 
ia UR ec =z D. 
7A dn, ur Sno 13 99a dt ( ) 


we see that the equation (A5.) still holds with the present more general marks of position of a 
moving system, and gives still the expression (B5.), supposing only, as before, that the two parts 
of the whole characteristic function are chosen so as to vanish with the time. 


It may not at first sight appear, that this rigorous transformation (B5.), of the partial 
differential equation (F.), or of the analogous equation (T.) with coordinates not rectangular, 
is likely to assist much in discovering the form of the part V, of the characteristic function V, 
(the other part V, being supposed to have been previously assumed;) because it involves under 
the sign of integration, in the term T, , the partial differential coefficients of the sought part V,. 
But if we observe that these unknown coefficients enter only by their squares and products, we 
shall perceive that it offers a general method of improving an approximation in any problem 
of dynamics. For if the first part V, be an approximate value of the whole sought function V, 
the second part V; will be small, and the term T, will not only be also small, but will be in 
general of a higher order of smallness; we shall therefore in general improve an approximate 
value V, of the characteristic function V, by adding to it the definite integral, 


i ` 
Vis Í (2-1) dr; (KS) 


though this is not, like (B5.), a perfectly rigorous expression for the remaining part of the 
function. And in calculating this integral (K5.), for the improvement of an approximation V,, 
we may employ the following analogous approximations to the rigorous formule (D.) and (E.), 


| i TIT US r, t, CRM al: 

6d, Ev $a, MyM; eee Sd, n "n? 

eV, PE. A F eV, Pi 5 

Sb, mbi; 8b, — Maba; ... Fr —m, 6); (L 3. 

9V, , 9V, L eV, , 

€ — M,C}; $e — Mg Co; ... Sent —My Cn; 

and y 

V, AA 5 
sg x 


or with any other marks of final and initial position, (instead of rectangular coordinates,) the 
following approximate forms of the rigorous equations (S.), 
eV, 8%, 98V, oT, Ky. 8T, 
le. o Beg! Be BT = Baga’ 
together with the formula (M5.); by which new formule the manner of motion of the system 
is approximately though not rigorously expressed. 


(N5.) 
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It is easy to extend these remarks to problems of relative motion, and to show that in such 
problems we have the rigorous transformation 


t 
Yaf (T, — T, T 4)dt, (05.) 
0 
and the approximate expression 
t 
Va= | (7, -7,)dt (pe) 
0 


V, being any approximate value of the function V, of relative motion, and V,, being the 
correction of this value; and T',, T',, being homogeneous functions of the second dimension, 
composed of the partial differential coefficients of these two parts V,,,V,., in the same way as 
T, is composed of the coefficients of the whole function V, . These general remarks may usefully 


be illustrated by a particular but extensive application. 


Application of the foregoing method to the case of a Ternary or Multiple System, 
with any laws of attraction or repulsion, and with one predominant mass. 


20. The value (68.), for the relative living force 27’, of a system, reduces itself successively 
to the following parts, 27°, 27’, ... 277("-9, when we suppose that all the » — 1 first masses 
vanish, with the exception of each successively; namely, to the part 


ieu AL HR 24 C12), (132.) 
Mı + 
when only m,, Mp, do not vanish; the part 
mam 
(2... B (£2.42. [2 : 
2T: m, +m, (& + T6?) (133.) 
when all but ma, m,,, vanish; and so on, as far as the part 
2:50, race. 134. 
sant (Fabi tib Ga) (134.) 


which remains, when only the two last masses are retained. The sum of these n — 1 parts is not, 
in general, equal to the whole relative living force 27T, of the system, with all the n masses 
retained; but it differs little from that whole when the first » — 1 masses are small in comparison 
with the last mass m, ; for the rigorous value of this difference is, by (68.), and by ien ) (133.) 
(134.),* 


2T, — 2T —27®— ... — 2T 0-9 = a (rw P) 75 (T —T )4.. D un 


2-1 (T(n-0 — T ) 


tz E, -mmy (5i — Er)? + (91 — n)? + ea fj)": (135. 


an expression which is small of the second order when the n — 1 first masses are small of the first 
order. If, then, we denote by V (9, y 9, ... V (-9, the relative actions, or accumulated relative 


23m T,— 
Tha 


" | 25, Er m T0-E,m (È+ +Å?) = a E,m,m, [(£ £&*- ... toed 


2 . 
"e m. Xm, T0--2»,T- if (,m,) T,+27,-+ E, Mi My (és — ér)? + Nt + oh 
n My My 
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living forces, such as they would be in the n — 1 binary systems, (m, m,), (m, Mn), ... (m, 4m,), 
without the perturbations of the other small masses of the entire multiple system of n points; 
so that 


yo- Í '2TOd, VO= Í "27e dt, ... V 0-D= Í ‘opad, (Q5.) 
0 0 0 
the perturbations being neglected in calculating these » — 1 definite integrals; we shall have, as 


an approximate value for the whole relative action V, of the system, the sum V,, of its values 
for these separate binary systems, 


Vi2V(O94-VO9-...4- V0. (R5.) 
This sum, by our theory of binary systems, may be otherwise expressed as follows: 


mm,w® mm, w? Mp- m, UD 


1= HR (S5.) 
m, "m My Ma * Mn Th, i m My 
if we put for abridgement 
1) 
wD = 40 9 + Í ip dr o, 
roO 
ra) 
wD =hD 9 +. Í rodra, (T5.) 
ro 
saad 4 
w- zd h—0»0 g—1 TT Í y'(n-—1) dr”-D, 
rn) 
In this expression, 
Aja» 
rs a faim ign 
NL (U5.) 


Searels LINDE mid 
p) = ES y 2 (m, 4 + m,) f"—» + 29-5 P^ go»? ? 


1D, ... 1—D being abridged expressions for the distances r™, ... 7-7, and f®, ... f^-? being 
abridgements for the functions f(^9), ... f-+™, of these distances, of which the derivatives, 
according as they are negative or positive, express the laws of attraction or repulsion: we have 
also introduced 2n — 2 auxiliary quantities h®, g®, ... 4("—D, g”—, to be eliminated or determined 
by the following equations of condition:* 


rD $0) 
0-99. | s dr®, 


ro) SAY 


5 op 
0290. Í a S) dr», (V5.) 


ra- Sp/n—D 


0=90-D4 Í or det, 


rr SAOD 


* [The equations (V5.), (X5.) are obtained by equating to zero the partial derivatives of V with respect to the 
arbitrary constants h®, g®. For this method of obtaining the characteristic function see Appendix, Note 2, p. 613.] 
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y 0) dr» (2) dr® in—1) dro») 
m Í ro) 70 =f ro m =f (n1) 7! (n—2)? (W5.) 
Sw) Bw Bu/n-» 
or $0 ag $g9— w= oy " (X5.) 
along with this last condition, 
mg V mg? omg» o0 fag o H,. avi 


Mi +My Ma + My Mg + My M ma a d My 


and we have denoted by 9”, ... 9"-D, the angles which the final distances 1%, ... r™—», 

of the first n— 1 points from the last or nth point of the system, make respectively with 

the initial distances corresponding, namely, rf”, ... r{”-». The variation of the sum V, is, 

by (S5.), 

M Mp WD MyM, ow? Mp1 Mp 8u (D 
Mı +My Ma +My E: m, a m, 


ôV, = ; (Z5.) 
in which, by the equations of condition, we may treat all the auxiliary quantities h®, g®, ... 
h”), g(^-D as constant, if H, be considered as given: so that the part of this variation 8V,,, 
which depends on the variations of the final relative coordinates, may be put under the 
form, 


Bum, Sut, Du? 
Mamie, (36 By, ogg, 8) 
mm, [SWP Sw?) , Swi? )e- 
vx Lm S ae t èna UR ic (A*.) 
m, ,m, (Sw) sur) ) 
m, 4 +My, (Se eat ae 1 i9 DOLO. sga m 


By the equations (T*.), (U5.), or by the theory of binary systems, we have, rigorously, 


1X2 1)\ 2 1)\ 2 
Ca) Ce i meme 


G6.) * (Spe) + Go.) Teatme By 


Su» 2 Su —»2 Spln—D\ 2 
F Suec — es —————— = 2 (n—1) n-D. 
( 865a | M ( 99,1 | id ( 8E | (m, 4 + Mp) fC +.2¢ 


and the rigorous law of relative living force for the whole multiple system is 


T, =U+H,, (50.) 
in which 
U =m, (Myf + mf + ... Em, 4/07») +2, mm SfE”, (C*.) 
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1/1 1 eV Y? (8VA?  [8V y? 
- T, (mtm) Se.) 652) * (62) | 
+3 nm) oe) * (5x) * (52) > 
2\m, m,) (Nos 875 dle 
eT 
2\My1 m 98, 90,1 90, 1 
E zit eV, dV, dV, V, x) 


SE, 8, SygSnn 80,00; 


We have therefore, by changing in this last expression the coefficients of the characteristic 
function V, to those of its first part V,,, and by attending to the foregoing equations, 

ee) as) 
M,+m,m,+m,\ SÉ; 96, Sq 0m, 56; 90,7 ; 


. (DS). 


T,2m,X,.m;f9-X-H,-t-m,X,. 
and consequently 


T T jy sui [re P Mn (s Sw Sw Sw ww (9? =) (Fe ) 
AG ORT 4k ET s 


(m, +m) (Ma + my) \ Sé; 8£, — 99; Sy, — 90; 9%; /) 


The general transformation of the foregoing number gives therefore, rigorously, for the 
remaining part V,, of the characteristic function V, of relative motion of the multiple system, 


the equation 
t Sw Sw = Sw Sw sw Sw 


t ; 
Y.-[ T Gdt4- X, mim; MES SE: one oe a dt; (G*.) 
0 
0 


E 
ui (m, +M) (My, + My) 
n 


EU and hence ü-l hy ly. 


M; d&; m, ÒE 
T,24X, — Lay gk POE ete.] 


t [T4 is formed from (D*.) by replacing V, by V,,. For the binary system composed of the points m;, m, we 
have 


* [From (X!.), 


expression (69) for T', gives 


Consequently the 


(OQ Mn (ere tay pre Nm ats OY 
2T; ws tax Vi ET Bde 


Ovi _ mmu, & UB. d 
és mtm, ee eee 


i i i i) 2 ` 
. op Tct m, (C (Ey oV; y) ELS e y (2 KER) } 22mm, 0. 2H19, 
"M mmy 0&; n Oni H ( a; m,4-m, \\ ôE: 3 Oy; ) 26 f 


mi--m 
where ————^ Hi=g, 
T4 T 


From equations (X!.), p. 126, 


\ 


Tjí-iL 


Mi +My oven? avee sav Se EE LA EM 
Minds ~ Ca) * Cx.) (CR) bi ma^ P 0f, ^ Om Om, li Obr 


ne (022 VR MT T 


= ar) + E n mh OO J Ot nO AE AR 
TE ry t\ On) tho BC ama Au On On. OG Oy 
=m,&,m; f+H,+... T,-Xm,m, fi +H,, ete. ] 


20 
HMPII 


www.rcin.org.pl 


154 IL. GENERAL METHOD IN DYNAMICS [20 
and, approximately, the expression 
j 1 "y 
V, 7X, mm, fi Spe inni + D |t: (Ht, 
n 


with which last expression we may combine the following approximate formule belonging in 
rigour to binary systems only, 
ów? , 8, 5w? 


(===, eL Es. I$. 
& sé, Ni 85; 6 8t, ( ) 
r2 Sw hae Sw ES Sw i 
Mise H , cai n= E. (K*.) 
| buo 
and = 590" (LS.) 


We have also, rigorously, for binary systems, the following differential equations of motion 
of the second order, 


y fO : fO 

= (my +m) a qom emos C= (my +m) or (at. 
96; oL; 

which enable us to transform in various ways the cana expression (H*.). Thus, in the 


case of a ternary system, with any laws of attraction or repulsion, but with one predominant 
mass ma, the disturbing part V,» of the characteristic function V, of relative motion, may be put 


under the form V,—m,m,W, (N®.) 
in which the coefficient W may approximately be expressed as follows: 
; 1 
wf {pom Ernie o a (0*) 
sf)  $B&f0 , fm 
orthus* — W-[ " er tm agr) 
dris sw) Su Sw) wD Bw? qu 
xl Mam WE t Th 8t, TENA (sg tns: 
(2) (2) (2) ` 
or finally, W= =Í ( nm TA in r ni a òf PEA be afe IL 
vB, (Q*.) 
"S. ów? . 5w9 SWD Sw 8w? j 
estes di^ hie: 7e 338 t$) 
In general, for a multiple system, we may put 
and approximately, V,,— E, . mm, WO; (R^.) 
© À) (i) 
Wi) — d fece de e “+m e+ A of AL 
: 6 
£ ers sw at sw sw dw “| P 
xl ksr Sé; Tt — 85, ksr. 8t, tue, Be 5B; TtYk y, , 
(k) (k) k) 
T anf (; ME 4 +m + +p AK 
(1*.) 


Th Sw hs Bu) Bu Saw 
A atait the tH Ba TP BB +H ya) 


* [Integrating by parts.] 
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Rigorous transition from the theory of Binary to that of Multiple Systems, by means 
of the disturbing part of the whole Characteristic Function; and approximate 
expressions for the perturbations. 


21. The three equations (K9.) when the auxiliary constant g is eliminated by the formula 
(L°.) are rigorously (by our theory) the three final integrals of the three known equations of the 
second order (M°.), for the relative motion of the binary system (m; m,,); and give, for such a 
system, the three varying relative coordinates £;, 5;, ¢;, as functions of their initial values 
and initial rates of increase «;, Bi, Yis «;. Pi» y;, and of the time ¢. In like manner the three 
equations (I9.), when g is eliminated by (L*.), are rigorously the three intermediate integrals of 
the same known differential equations of motion of the same binary system. These integrals, 
however, cease to be rigorous when we introduce the perturbations of the relative motion of 
this partial or binary system (m,;m,,), arising from the attractions or repulsions of the other 
points my, of the whole proposed multiple system; but they may be corrected and rendered 
rigorous by employing the remaining part V,, of the whole characteristic function of relative 
motion V, , along with the principal part of approximate value V , . 


The equations (X!.), (Y!.) of the twelfth number give rigorously 
DOR UI 1 8V, 1, V, i 1 oY, TAG, 


= d APIS sieut AME r sai? AEF eee, 6 
Mu E ee A gear mar, (UM 
and 
: RIAM B nor, p. Ele. Aon NE, KXIW ho OF, : 
TUM m; "Head ni “P= t a m," 8B; Ym, Sy; ea P (V®.) 
and therefore, by (A°.), 
sw) my owe 1 OV. 1 39V, 
AE nH 
9€; My +My dg; Mi dg; My o£; 
Sw m, sw 18, 1.88 
=7;,—,- ET MUSEUM M el Ws. 
Sy My +My Syn, Mi Sy; m, | dy; uk. 
bw V oy Me Sw 13V, 1 BV, 
ob; om, tm, OF, mob m, '8Lb 
and similarly 
Sw _ yoke a Sw) JY. 1 V, 
a, 5 amy tm, Baty mde, m, "Bay 
sw) m, sw 18, QA 8V,, > 
CM CP eru em, SBa Mm, BB, m, BB,’ dx 
ówO' , m7 uM" a 1. SF, 
by, oe ` My + My Sy, tm E n, ^ y; ' 


the sign of summation £, referring only to the disturbing masses m; , to the exclusion of m; and 
m, ; and these equations (W?.), (X9.) are the rigorous formule, corresponding to the approximate 
relations (I9.), (K9.). In like manner, the formula (L.) for the time of motion in a binary system, 
which is only an approximation when the system is considered as multiple, may be rigorously 


20-2 
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corrected for perturbation by adding to it an analogous term deduced from the disturbing part 
V, of the whole characteristic function; that is, by changing it to the following: 


Su BV, 


Is ai $6) - 
(= 300 * 3H” (Y » 
which gives, for this other coefficient of w, the corrected and rigorous expression 
Su , BV. 6 
8g? -— t OF ôH, (Z .) 


V, being here supposed so chosen as to be rigorously the correction of V,,. If therefore by the 
theory of binary systems, or by eliminating g between the four equations (K*.) (L9.), we have 
deduced expressions for the three varying relative coordinates £;, 4;, 6; as functions of the 
time ¢, and of the six initial quantities «;, B;, y;, «;, Pi» Yis Which may be thus denoted, 


éi = ha (e. Bis yi oi Pio Yi» t)» 
Ni = $a (e; Bis Yis tis Bis Yi» t)» (A7.) 
Ci = Ps (5 Pis Yi» ti Bio vist) i 
we shall know that the following relations are rigorously and identically true,* 
TA: Yi bois hero ee eus] 
Ye Ba? SB B MOD 
AR e AMA den A) 
Ni — P2 | Xis Pi» Yi» Sa,” 5B,” y," dg)” 
Gods Bev ME. asi Mam iuo]. 
ba,’ ~ 8B,’ Dy, yo 


(B*.) 


and consequently that these relations will still be rigorously true when we substitute for the 
four coefficients of w their rigorous values (X*.) and (Z?.) for the case of a multiple system. We 
may thus retain in rigour for any multiple system the final integrals (A’.) of the motion of a 
binary system, if only we add to the initial components «;, B;, y; of relative velocity, and to 
the time ż, the following perturbational terms: 


my SW, 1 Wg, 1 Va 


cme dT omo da, me ' 8a,” 
Afi eX, E at = Ma ent si (C7) 
) m ou E E Vs V. 
ix Avi 3, "my em, Sye Bye. mi "d a [x 
At= -3H (D'.) 


* [Otherwise we would get one or more equations connecting the initial and final points. Of course when 
the rigorous values are substituted, £;, ;, £; denote the actual coordinates at the time ż.] 
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In the same way, if the theory of binary systems, or the elimination of g between the four 
equations (I*.) (L°.), has given three intermediate integrals, of the form 
& = pa (Éis Mi os Bio Yi» t)» 
Ni = Pa (Ei 05 Cis tis Bis Yis t), (E?.) 
C = Pa (Eis Mis Cio tis Bio Yis t), 
we can conclude that the following equations are rigorous and identical, 


i) 
a Py (ton toa Boves. 
i 
Se= Pe (Eom ne Boves]. vl 


c, mh (Feo te tis Bi Yi» asl 


and must therefore be still true, when, in passing to a multiple system, we change the coefficients 
of w to their rigorous values (W*.) (Z*.). The three intermediate integrals (E?.) of the motion 
of a binary system may therefore be adapted rigorously to the case of a multiple system, by 
first adding to the time ¢ the perturbational term (D’.), and afterwards adding to the resulting 
values of the final components of relative velocity the terms 


m, Ow? 18V, SV o 
Agi Ug, v Tc TE 
5-Y, my + m B. Um 9€, : a 
MERECE T y an 


Sd ee Ce CFR" MD 
mM,+mM, Sn, món; m, ' SH; 
my 3B. 1 Bia’, 1. BV 
““my +m, 95, m; 5; m, ' Òb 
22. To derive now, from these rigorous results, some useful approximate expressions, we 
shall neglect, in the perturbations, the terms which are of the second order, with respect to the 
small masses of the system, and with respect to the constant 2H, of relative living force, which 
is easily seen to be small of the same order as the masses: and then the perturbations of the 
coordinates, deduced by the method that has been explained, become 
i ò$; O84 a.» , O88 
Ag, = Zon METÀ ABi+ 5 Avi t gg At 
i Oni A gr OM anr LOI 
cs Aa m ABL so Ay At, H’. 
Sa, i 8 Bi B; Sy} Yi ot ( ) 
4 8L; 8L, at 
AC, = 5. Aai + zs ABL sos Ay + At, 
C= Bar ag; Pct; An 8t 
in which we may employ, instead of the rigorous values (C’.) for Ax, AB;, Ay;, the following 
approximate values: 


Ati== 


An= 


K 


AB, EX, — M—— te (I*.) 
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To calculate the four coefficients 

Vig Vig Vig BV, 

da,’ dB,’ dy,’ dH,’ 
which enter into the values (I’.) (D7.), we may consider V,,, by (R*.) (T9.), and by the theory of 
binary systems, as a function of the initial and final relative coordinates, and initial components 
of relative velocities, involving also expressly the time ¢, and the » — 2 auxiliary quantities of 
the form g?; and then we are to consider those initial components and auxiliary quantities and 
the time, as depending themselves on the initial and final coordinates, and on Z,. But it is not 
difficult to prove, by the foregoing principles, that when ¢ and gy are thus considered, their 
variations are, in the present order of approximation,* 


zo) eei a 
REES Aa SAG, PC (K?.) 
EE : 
2,.m gà 
jiu fea Bul) 
ad bg = (13) (3-8, ws]: (L7.) 


the sign of variation 9, referring only to the initial and final coordinates; and also that 
Sw SE, Fw SSi Pw Sé,  d'w Sé, 
39 dt Bx, 8g Bai * 55,7088, * 5y, 8g Sy, 
along with two other analogous relations between the coefficients of the two other coordinates 
Ni» Či; from which it follows that t and gy, and therefore a; , B; , y; , may be treated as constant, in 
taking the variation of the disturbing part V,,, for the purpose of calculating the perturbations 
(H*.): and that the terms involving A£ are destroyed by other terms. We may therefore put 


(M’.) 


simply} Aé,= oe déi A Bi + d I " 

1 òb; i , 

s E 
Aq, sl d ABi+ E Ayi» (N’.) 

à 

Al uis a ABA 

* [We have iust But w™® depends on the initial and final coordinates and also on g™, Thus 
Sw By OD 


ago sg Tho 


(L’.) follows immediately and then (K’.) by summation. Equation (M".) is obtained by differentiating (B’.) partially 
with respect to g™.] 
T [To obtain (N’.): Vj, — Z,m,m, WC 9, where WG, ®© is given by (T*.). Three groups of terras come from 


) 
(i) From the integrated part of W^? we get "Ws v - This cancels against the first term of Ao; in (I’.). 
k 


(ii) From ¢ which occurs “expressly” (explicitly) in the upper "T of the integral. From (K’.) 
E: Sw Sw Sw 8*5) 71 
Bar magos) saci / = (Gu) "al m= In (5) 


Since At= — SH , we see, by (M’.), that this group of terms is cancelled by the terms involving At. 
(iii) The remaining terms which appear in (O".).] 
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employing for A«; the following new expression, 
M IERCH Sah (18 RG | 
Aa; —XÀ, mf 5 ia ot owe dt 
ME § RG K) à Bi SRG 
ba; Jo dB; Ba; Jo By 
together with analogous expressions for Af;, Ay;, in which the sign of summation X, refers to 
the disturbing masses, and in which the quantity 
sf — af af% 
XA Ni 85, b 8C, 
is considered as depending on «;, Bis Yis «;; Bis Yis Aro Bas Yro Ako Bas Yh» t, by the theory of binary 
systems, while «;, B;, y; are considered as depending, by the same rules, on «;, Bis Yio Eis Mis Gi 
and f. 


It may also be easily shown, that* 


Mei, MM dEr Sot M e. 
Òa; 5a; 9B; 9B; oy; Oy; 8a; 


with other analogous equations: the perturbation of the coordinates £; may therefore be thus 


(07.) 


RiD=f6V4E LET) 


expressed, 
T o£, [ORGO Sé, (t5. R00 
Ae E, ms. óc; a- s |, à; a 
8£; [ED £, — 
<s =s t- | =r l R”. 
5B; Jo 8B, “73B: o 3P; jan 
E Rr Ec MES | 
óy;Jo Oy; ôy; Jo dy; i 


and the perturbations of the two other coordinates may be expressed in an analogous manner., 


Tt results from the same principles, that in taking the first differentials of these perturbations 
(R’.), the integrals may be treated as constant; and therefore that we may either represent the 
change of place of the disturbed point m;, in its relative orbit about m, , by altering a little the 
initial components of velocity without altering the initial position, and then employing the 
rules of binary systems; or calculate at once the perturbations of place and of velocity, by 
employing the same rules, and altering at once the initial position and initial velocity. If we 
adopt the former of these two methods, we are to employ the expressions (O".), which may be 
thus summed up, 


9 ' RH dt 
| Se í " 
t 
Af -X, mag [Boras (S?.) 
8 (tv. 
Ayi—-X, m. | RG» dt; 


0a; Pw) opi Ox, Pwd əy 


OB, 00408, Om Oy, Oyða Üo” 


* 
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and if we adopt the latter method, we are to make, 


t i, k) t i, k) 
| ie dt, As,- -X,-m, | Mr d 
0 0 
t § RG t § RO. 
-—2X,. ——,— dt, T 
Aem de mJ. pp, edad ETE mf 5B; diii 


t6 RG) t § RG) 
Ay; = i m, |. Sy; dt, Ay; = —X, sm, | hz 


The latter was the method of Lagrange: the former is suggested more immediately by the 
principles of the present essay.* 


General introduction of the Time into the expression of the Characteristic Function 
in any dynamical problem. 
23. Before we conclude this sketch of our general method in dynamics, it will be proper to 
notice briefly a transformation of the characteristic function, which may be used in all applica- 
tions. This transformation consists in putting, generally, 


V=tH +8, . (U*.) 
and considering the part S, namely, the definite integral 
t 
S= Í (T+U)dt, (V7.) 
0 


as a function of the initial and final coordinates and of the time, of which the variation is, by our 
law of varying action, 


8S — — H8t -- Xi. m (x'8x —a'8a +y’ dy — b'8b 4- 2/82 — c'8c). (W’.) 
The partial differential coefficients of the first order of this auxiliary] function S are hence, 
òS | 
mandi: 7 
x7 0H (X*.) 


oS * 
By, is m, 450 (Y’.) 


* [If we give variations A«;, Af,, gU the new variations of «;, ;, y; will be 


ba] Saf 
Ace, + 2% sei An + p MB, te Ma Fe An et AB Se ^u 
since 9B; uL dw 
8x, 58; Sopp 


where Aa; etc. mean the values given in (07.). By giving the values of Ax;, A8;, Ay; in (T’.) we get the values of 
Aa’, in the same set of equations. See Appendix, Note 3, p. 622. Cf. Lagrange, Mécanique Analytique, Tome 1, 
2nd Part, Section v. Lagrange, “Sur la théorie générale de la variation des constantes arbitraires," Mémoires de 
l'Institut (1808), p. 257. For a comparison of the methods of Lagrange, Poisson and Hamilton in the Theory of 
Perturbations, see Lovett, Quarterly Journal of Mathematics, xxx (1899), pp. 47-149.] 

T [Hamilton perceived at this stage that it was inconvenient to have H appearing in the characteristic function; 
so he transformed to another function which contained ¢ as a variable instead of H. The transformation to adopt is 
obvious from the following equations: 

8 V =Em (a'0x — a'8a + y/8y — b'8b -- 2/82 — c/8c) -t8H, 
~ 6(V — Ht) Xm (z'dx —a'da 4 ...) — Hàt.] 
1 [Hamilton changes the name “auxiliary” to “principal” in the Second Essay.] 
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8S ^ Re "ug 3 : $ 
and ET S 5B, 0e E^ — ne; . (Z?.) 
These last expressions (Z’.) are forms for the final integrals of motion of any system, corre- 
sponding to the result of elimination of H between the equations (D.) and (E.); and the expres- 
sions (Y’.) are forms for the intermediate integrals, more convenient in many respects than the 


forms already employed. 


24. The limits of the present essay do not permit us here to develope the consequences of 
these new expressions. We can only observe, that the auxiliary function S must satisfy the 
two following equations, in partial differentials of the first order, analogous to, and deduced 
from, the equations (F.) and (G.): 


ORORO a 


Za (n ^ 


and that to correct an approximate value S, of S, in the integration of these equations, or to 
find the remaining part $,, if 


and 


S - 8,-- S,, (C*.) 
we may employ the symbolic equation 
d 8 1/888 888 888 
dà Lira bp Be fi 
which gives, rigorously, 
ee ac m 
if we establish by analogy the definition 
olen slay ME) m 
and therefore approximately 
saf: (U—U,) dt, (G$.) 


the parts S,, S, being chosen so as to vanish with the time. These remarks may all be extended 
easily, so as to embrace relative and polar coordinates, and other marks of position, and offer 
a new and better way of investigating the orbits and perturbations of a system, by a new and 
better form of the function and method of this Essay. 


\ 


March 29, 1834. 
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